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B : ABSTRACT 

We investigate the eilects of thermonuclear reaction rate variations on ^^Al production in 
' massive stars. The dominant production sites in such events were recently investigated by using 

stellar model calculations: explosive neon-carbon burning, convective shell carbon burning, and 
convective core hydrogen burning. Post-processing nucleosynthesis calculations are performed 
Pm I for each of these sites by adopting temperature-density-time profiles from recent stellar evolution 

C/^ ■ models. For each profile, we individually multiplied the rates of all relevant reactions by factors of 

(~| I 10, 2, 0.5 and 0.1, and analyzed the resulting abundance changes of ^^Al. In total, we performed 

C^H' ~ 900 nuclear reaction network calculations. Our simulations are based on a next-generation 

nuclear physics library, called STARLIB, which contains a recent evaluation of Monte Carlo 
\—( • reaction rates. Particular attention is paid to quantifying the rate uncertainties of those reactions 

c/3 I that most sensitively influence ^^Al production. For stellar modelers our results indicate to what 

degree predictions of ■^^Al nucleosynthesis depend on currently uncertain nuclear physics input, 
while for nuclear experimentalists our results represent a guide for future measurements. We also 
investigate equilibration effects of ^^Al. In all previous massive star investigations, either a single 
I species or two species of ^^Al were taken into account, depending on whether thermal equilibrium 

was achieved or not. These are two extreme assumptions and in a hot stellar plasma the ground 
and isomeric state may communicate via 7-ray transitions involving higher-lying ^^Al levels, 
ly-^ ■ We tabulate the results of our reaction rate sensitivity study for each of the three distinct mas- 

sive star sites referred to above. It is found that several current reaction rate uncertainties influ- 
ence the production of ^®A1. Particularly important reactions are ^^Al(n,p)^^Mg, ^^Mg(a,n)^*Si, 
^^Mg(n,7)^^Mg and ^•^Na(a,p)^^Mg. These reactions should be prime targets for future mea- 
surements. Overall, we estimate that the nuclear physics uncertainty of the ^^Al yield predicted 
by the massive star models explored here amounts to about a factor of 3. We also flnd that 
taking the equilibration of ^^Al levels explicitly into account in any of the massive star sites 
, ^ investigated here has only minor effects on the predicted ^^Al yields. Furthermore, we provide 

■ for the interested reader detailed comments regarding the current status of certain reactions, 

including i2c(i2c,n)23Mg, ^^Na{a,p)^^Mg, ^^Mg{a,n)^^Si, 26Al'"(p,7)27si, 26Al(n,p)26Mg and 
26Al(n,a)23Na. 

Subject headings: gamma rays: theory — nuclear reactions, nucleosynthesis, abundances — stars: evo- 
lution — supernovae: general 
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1. Introduction 

The radioisotope ^^Al is of outstanding impor- 
tance for 7-ray astronomy and cosmochemistry. It 
has been discovered in three distinct sites: (i) in 
the Galactic interstellar medium via detection of 
its decay emission line at 1809 kcV (Mahoncy ct 
al. 1982, Diehl et al. 1995); (ii) in meteorites via 
observed excesses of its radioactive decay (daugh- 
ter) product ^^Mg (MacPherson, Davies & Zin- 
ner 1995), implying an injection of live ^^Al into 
the early Solar System nebula; and (iii) in preso- 
lar dust grains, again via detected ^^Mg excesses 
(Hoppe et al. 1994; Huss, Hutcheon & Wasserburg 
1997), that are uncontaminated by solar system 
material and thus are of likely stellar origin. Iden- 
tification of the main sources of ^^Al would have 
far-rcac;liing implications, ranging from questions 
related to the circumstances and conditions of the 
Solar System birth to imposing strong constraints 
on the chemical evolution of the Galaxy. A num- 
ber of different sources have been suggested over 
the years: AGB stars, classical novae, Wolf-Rayet 
stars, and core collapse supernovae. For reviews, 
see Prantzos & Diehl (1996) or Diehl & Timmes 
(1998). However, the origin of ^^Al remains con- 
troversial. 

The observation of Galactic 7-rays from ^^Al 
is important since it provides unambiguous direct 
evidence for the theory of nucleosynthesis in stars. 
The half-life of ^^^Al amounts to 7.17 x 10^ y and 
is small compared to the time scale of Galactic 
chemical evolution (f»10^° y). Consequently, nu- 
cleosynthesis is currently occurring in the inter- 
stellar medium and, in particular, ^^Al is synthe- 
sized throughout the Galaxy. Prom the observed 
7-ray intensity, depending on the assumption for 
the density distribution, a present-day ^®A1 equi- 
librium mass of « 2 — 3 Mq in the entire Galaxy 
has been inferred (Diehl et al. 2006). The obser- 
vational evidence favors in this case massive stars 
as a source: first, the all-sky map of the 1809 
keV 7-ray line detected by the COMPTEL instru- 
ment onboard CGRO showed that ^^Al is confined 
along the Galactic disk and that the measured 
intensity is clumpy and asymmetric (Pliischke et 
al. 2001); second, the comparison between the 
26 Al all-sky map from COMPTEL to other all-sky 
maps for different wavelengths (Knodlseder 1999) 
revealed that the 1.8 MeV 7-ray emission is cor- 



related with the Galactic free-free emission, which 
traces the distribution of ionized gas in the inter- 
stellar medium, observed in the microwave domain 
by the COBE satellite; third, the measurement of 
the 1809 keV line Doppler shift by the SPI spec- 
trometer onboard INTEGRAL demonstrated that 
26 Al co-rotates with the Galaxy and hence sup- 
ports a Galaxy- wide origin (Diehl et al. 2006). 
Here we investigate the bulk production of ^^Al 
in massive stars. A study of ^^Al/^'^Al ratios ob- 
served in meteorites and presolar grains will be 
subject of a separate study. 

Massive stars may produce ^^Al during sev- 
eral different phases of their evolution: (i) dur- 
ing pre-supernova stages in the C/Ne convective 
shell, where a fraction of the ^^Al survives the 
subsequent explosion and is ejected into the inter- 
stellar medium (Arnett & Wefel 1978); (ii) during 
core collapse via explosive Ne/C burning (Arnett 
1977), where the ejected ^^Al yield may perhaps 
be modified by the i/-process via neutrino spal- 
lation (Woosley et al. 1990); and (iii) in Wolf- 
Rayet stars, i.e., stars with masses in excess of 
about 30 Mq, which experience such a strong mass 
loss that even layers located within the H con- 
vective core, hence significantly enriched in ^^Al, 
are ejected into the interstellar medium (Pala- 
cios et al. 2005). These ^^Al production mech- 
anisms (with the exception of the z/-process) were 
recently analyzed in detail by Limongi & Chieffi 
(2006) by using extensive stellar model calcula- 
tions of solar metallicity stars in the mass range 
of UMq < M < 12OM0. In that work they also 
emphasized the impact of rate uncertainties for 
selected reactions on the final ^^Al yields. An- 
other discussion along these lines can be found in 
Woosley & Heger (2007). 

In the present work we expand this effort by 
presenting a comprehensive investigation of the 
impact of nuclear reaction rate uncertainties on 
the synthesis of ^^Al. Our method is similar in 
spirit to earlier nuclear reaction rate sensitivity 
studies that addressed the nucleosynthesis in clas- 
sical novae (Iliadis ct al. 2002) and type I x-ray 
bursts (Parikh et al. 2008). The general strat- 
egy consists of varying the rates of many reactions 
by different factors (in this work, 10, 2, 0.5 and 
0.1) and to analyze the impact of each individ- 
ual reaction rate change on the final ^^Al yields. 
Once the yield changes are established for this 
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grid of rate variation factors, more realistic abun- 
dance changes based on actual rate uncertainties, 
if available, are considered. At present it is not 
feasible to perform this computationally inten- 
sive procedure with a self-consistent stellar model. 
Instead, we extract representative temperature- 
density-time profiles from recent stellar evolution 
models of massive stars and execute a large num- 
ber of post-processing reaction network sensitivity 
calculations using these profiles. Our goal is two- 
fold. On the one hand, we would like to quantify to 
what degree predictions of ^^Al yields depend on 
currently uncertain nuclear physics input. On the 
other hand, by identifying the "most important 
nuclear reactions", our results represent a guide 
for future measurements. 

There are a number of novel aspects about the 
present work. First, we employ a new- generation 
library of nuclear reaction and weak interaction 
rates, called STARLIB. This library is partially 
based on a recent evaluation of experimental 
Monte Carlo reaction rates (Iliadis et al. 2010). 
Besides recommended reaction rates for a grid 
of temperature values between 1 MK and 10 GK, 
the library includes in addition for many reactions 
the rate uncertainty factor at each temperature. 
In fact, this work represents the first application 
of STARLIB. Second, we carefully investigate the 
equilibration effects of Al. At least two species of 
^^Al take part in the nucleosynthesis, the ground 
state and the isomeric state. In all previous in- 
vestigations, either a single species or two species 
of ^^Al were taken into account, depending on 
whether thermal equilibrium is achieved or not. 
Obviously, these are two extreme assumptions and 
in a hot stellar plasma the ground and isomeric 
state may "communicate" via 7-ray transitions 
involving higher-lying ^^Al levels. 

Our study has some obvious limitations. First, 
since we perform post-processing calculations, we 
necessarily focus our investigation on the effects 
of nuclear reaction rates. In other words, the im- 
portant effects of convectior0, mass loss, rotation. 



^It can be shown analytically that the abundance evolution 
in a convective region in which the turnover time is fast 
enough to ensure a flat abundance profile of the various 
nuclear species is equivalent to the evolution of a single 
mesh in which each local thermonuclear rate, (cri;) is re- 
placed by its mass-weighted average over the convective 
region, i.e., Y2k ('^^)ijk dmk/mtot- Hence it is perfectly 



and so on, are outside the scope of the present 
work. This also implies that our simulations are 
unsuitable for defining absolute ^^Al yields. In- 
stead, we claim that our procedure is useful for 
exploring the effects of ^^Al abundance changes 
that result from reaction rate variations. Second, 
we only explore a few temperature-density-time 
evolutions that are representative of solar metal- 
licity stars. A more comprehensive study covering 
a broad range of stellar masses and metallicities 
is also beyond the scope of this work. Third, it 
is well-known that the radioisotope ^"Fe (half-life 
of 2.62 X 10^ y) is likely co-produced with ^^Al in 
massive stars and that their abundance ratio pro- 
vides a sensitive constraint on stellar models (see 
Limongi & Chiefh 2006; Woosley & Heger 2007; 
and references therein). Indeed, Galactic 7-rays 
from the decay of ^°Fe have been detected by both 
RHESSI and the SPI spectrometer onboard INTE- 
GRAL, and the observed 7-ray line flux ratio for 
'^"Fe and ^'^Al amounts to w 0.1-0.2 (Harris et al. 
2005). In the present work we only focus on the 
nucleosynthesis of ^^Al and leave a similar sensi- 
tivity study for ^°Fe to future wor Finally, by 
individually varying each rate and leaving all other 
rates at their nominal values, we disregard any cor- 
relations among different reactions. In our opin- 
ion, no single study will cover all of the possible 
uncertainties, but each approach has advantages 
and disadvantages. We present only one realiza- 
tion of a sensitivity study, similar to the proce- 



plausible to use a single point evolution for a convective re- 
gion if the turnover time is faster that the nuclear burning 
time (e.g., for core H burning). Furthermore, we performed 
the calculation using un-weighted rates instead of mass- 
weighted average rates because we are mostly interested in 
^®A1 abundance changes. 
^In the recent work of Tur, Heger & Austin (2010), the im- 
pact of triple-o and ^■^C(o,7)^^0 rate variations on the 
^®A1, ^^Ti and ^''Fe yields were investigated using stellar 
evolution and explosion models. The authors claimed that 
"...over a range of twice the experimental uncertainty, a, 
for each helium-burning rate, the production of Al, ^"Fe, 
and their ratio vary by factors of 5 or more...". By looking 
in detail at their results, it is clear that the effects on the 
^''Fe yield are indeed large. However, it becomes also ap- 
parent that the effects on the ^®A1 yield are much smaller. 
For example, their 25Mq model and adopting the initial 
abundances from Lodders (2003) provides a factor of 1.5 
change in ^''Al yield if the triple-o and ^■^C(a,7)^®0 rates 
are individually varied by their Icr experimental uncertain- 
ties (see their Tab. 3). As will be seen, the rate uncertainty 
effects explored in the present work result in significantly 
larger ^^Al yield variations. 
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dure applied in Iliadis ct al. (2002) and Parikh et 
al. (2008). Interestingly, the latter work explored 
two methods: the one applied here and also a 
Monte Carlo procedure. It was found by Parikh et 
al. (2008) that very similar results were obtained 
with "...minor differences attributed to such corre- 
lation effects...". We feel that a Monte Carlo pro- 
cedure makes most sense if it is performed with 
reliable experimentally based reaction rate proba- 
bility densities. However, as will become apparent 
below, we do not have this information for all of 
the important reaction rates yet and thus leave 
such a study to future work when an update of 
STARLIB becomes available. 

This paper is organized as follows. Our proce- 
dure is explained in more detail in § 2, including 
a discussion of stellar models, the equilibration of 
^^Al, and a description of the library STARLIB. 
The results of reaction rate sensitivity studies are 
presented in § 3 for the three predicted main sites 
of ^^Al synthesis: explosive Ne/C burning, con- 
vective C/Ne shell burning, and convective H core 
burning. A simimary and conclusions are given in 
§ 4. More information on reaction and decay rates, 
together with a discussion of individual reactions, 
is provided in the Appendix. 

2. General Procedure 
2.1. Massive star models 

The stellar models adopted in the present work 

are those presented in Limongi & Chiefff (2006). 
For the sake of completeness, we summarize the 
basic properties of these models and the main re- 
sults concerning the production of ^^Al in mas- 
sive stars. The evolution of each stellar model was 
computed from the pre-main sequence phase up to 
the onset of the iron core collapse by using the stel- 
lar evolutionary code FRANEC (Frascati RAph- 
son Newton Evolutionary Code, release 5.050218). 
The kernel of this code has been presented in 
Limongi & Chiefff (2003) (and references therein). 
Here we will only mention recent updates. First, 
the convective mixing and the nuclear burning 
were coupled together providing a set of diffusion 
equations that are linearized and solved simulta- 
neously by means of a Newton-Raphson method. 
This coupling is extremely important in all sit- 
uations where the nuclear burning timescale of a 
given nuclide is comparable to the mixing turnover 



time. Thus the interaction between the local nu- 
clear burning and the convective mixing cannot be 
disregarded. The nuclear network adopted was the 
same as in Limongi & Chieffi (2003) and the ther- 
monuclear reaction rates were up-to-date at the 
time when these models were computed (see Tab. 
1 of Limongi & Chieffi 2006). The nucHde ^SAl was 
treated in a distinct manner, by assuming two sep- 
arate species (for the ground and isomeric state) 
for temperatures below T « 1 GK, and a sin- 
gle (thermalized) species above this temperature. 
Mass loss was included following the prescriptions 
of Vink et al. (2000) for the blue supergiant phase 
(Te// > 12000 K), de Jager et al. (1988) for the 
red supergiant phase {T^ff < 12000 K), and Nugis 
& Lamers (2000) for the Wolf-Rayct phase. All of 
these solar metallicity models had an initial He 
mass fraction of 0.285 and a global metallicity (by 
mass) oi Z = 0.02. The relative abundances for 
the various nuclear species were adopted from An- 
ders & Grevesse (1989). 

The explosion of the mantle of the star was 
started artificially, by instantaneously imparting 
an initial velocity of vq to a mass coordinate cor- 
responding to w 1 Air:) of the presupernova model. 
Such a mass coordinate relates to a region located 
well within the iron core and is chosen in such a 
way that the initial conditions should not affect 
the properties of the shock wave too much at a 
time when it approaches the Fe-Si interface. The 
formation and propagation of the shock wave, gen- 
erated in such a way, was calculated by means of 
a computer code that solves the fully compress- 
ible reactive hydrodynamic equations by applying 
the Piecewise Parabohc Method (PPM) of Colella 
& Woodward (1984), using a Lagrangian scheme. 
The chemical evolution of the matter was com- 
puted by coupling the same nuclear reaction net- 
work that was adopted in the hydrostatic calcu- 
lations to the system of hydrodynamic equations. 
The free parameter vq was properly adjusted in 
order to obtain a given final kinetic energy of the 
ejecta or, equivalcntly, to eject a given amount of 
mass above the Fe core. Since -^^Al was synthe- 
sized in regions relatively far away from the Fe 
core (see below), its final yield did not depend on 
the particular choice of vq provided that at least 
a minimum amount of ^^Ni was ejected. 

Based on this scit of presupernova models and 
simulated explosions, it was found that ^^Al is 
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produced in massive stars in three distinct evo- 
lutionary stages: core H burning, C convective 
shell burning just prior to the core collapse, and 
explosive Ne/C burning. Any ^^Al produced by 
these massive stars is ejected into the interstel- 
lar medium, both by stellar winds and the ex- 
plosion, in different proportions depending on the 
initial stellar mass. It was also found that ^^Al 
was mainly produced by explosive Ne/C burning 
over most of the initial mass interval between 11 
and 120 Mq. Only for the more massive stars, 
say, M > 6OM0, did the wind component (pro- 
duced during core H burning) become important 
(see Tab. 3 and Fig. 2 in Limongi & Chieffi 2006 
for details). In the present work, we based our 
study on 20, 60 and 80 Mq model stars for ex- 
ploring ^^Al yield sensitivities in explosive Ne/C 
burning, C convective shell burning and core H 
burning, respectively. These choices of model stars 
are motivated by the fact that they provide rela- 
tively large ^^Al yields. 

2.2. Nuclear physics library 

The nuclear physics input for the present post- 
processing studies is based on a new-generation 
library, called STARLIB. It originated from a pre- 
vious version of REACLIB (originally created by 
F.-K. Thielemann) that one of us modified over 
the years and was used for all of the reaction net- 
work calculations presented in Iliadis 2007. 

At that point in time several important changes 
occurred. A recent evaluation of reaction rates 
for the A=14-40 target range (Ihadis et al. 2010) 
was completed. These 62 experimental rates are 
based on a Monte Carlo technique, allowing for a 
rigorous definition of recommended reaction rates 
and their associated uncertainties. The Monte 
Carlo procedure also provides, for the first time, 
for any given temperature the (output) reaction 
rate probability density function that is based on 
the (input) probability densities of measured nu- 
clear physics quantities (such as S-factors, reso- 
nance energies, resonance strengths, upper limits 
in spectroscopic factors, etc.). From the cumula- 
tive distributions of the rate probability densities, 
a low rate, median rate and high rate can be de- 
fined as the 0.16, 0.50 and 0.84 quantiles, respec- 
tively, assuming a coverage probability of 68%. 
The meaning of these rates is in general differ- 
ent from the commonly reported, but statistically 



meaningless, literature expressions "lower limit", 
"nominal value" and "upper limit" of the total re- 
action rate. It is important to emphasize that the 
Monte Carlo rates incorporate both statistical and 
systematic uncertainties, as explained in detail in 
Longland et al. (2010). Furthermore, it has been 
shown in Longland et al. (2010) that in the major- 
ity of cases the Monte Carlo rate probability den- 
sity function can be approximated by a lognormal 
distribution, which is determined by only two pa- 
rameters: the lognormal location parameter /i and 
the lognormal spread parameter a. The former 
parameter determines the recommended reaction 
rate via Na (o'i')rec = while the latter param- 
eter corresponds to the rate factor uncertainty via 
f.u. = e''. 

The information on the rate probability density 
was not available previously and opens interesting 
windows of opportunity for Monte Carlo studies 
of nucleosynthesis and energy generation in stars. 
However, it becomes clear from the above dis- 
cussion that three quantities (T, Na {<^v)^^^, log- 
normal a) instead of the traditional two (T and 
Na (o'w)rec) need to be reported so that the user 
can calculate the rate probability density for each 
reaction at each temperature. Therefore, it was 
decided to convert our 2007 version of the REA- 
CLIB, which lists the recommended reaction rates 
as analytical functions of temperature by employ- 
ing a number of rate fitting parameters, to a tab- 
ular format. To be precise, the rate tables are 
directly derived from the fitting parameters and 
not from any tabular rates given in the original 
publications. The new format consists of three 
columns and lists for each reaction the tempera- 
ture, the recommended rate and the rate factor 
uncertainty on a grid of 60 temperatures between 
1 MK and 10 GK, allowing for an accurate in- 
terpolation between grid points. At this stage 
the rate factor uncertainty for each reaction is set 
equal to a nominal value of 10. In a subsequent 
step, the rates and factor uncertainties of 62 re- 
actions in the A= 14-40 region were replaced with 
their exact Monte Carlo results. In addition, the 
rates of the following interactions were replaced 
with more recent informatiorH: (i) 10 Big Bang 



^Assuming that the reaction rate probability density func- 
tion can be approximated by a lognormal distribution, it 
can be shown (Iliadis et al. 2010) that for a coverage prob- 
ability of 68% the lognormal spread parameter is given by 
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reactions, using the rates of Descouvemont et al. 
(2004), which were derived from an R-matrix de- 
scription of the available data; (ii) 30 reactions 
from the NACRE evaluation of experimental rates 
(Angulo et al. 1999), in the mass range of A=l- 
26; (iii) (n,7) reactions based on the KADoNiS 
vO.2 evaluation of experimental rates (Dillmann 
et al. 2006); (iv) a number of special reactions, 
such as "N(p,7)i5 (Bertone et al, m prepara- 
tion) and i2c(a,7)i60 (Kunz et al. 2002); (v) 
550 experimental rates for /3-decays and /3-delayed 
particle decays including associated uncertainties, 
calculated from the half lives and branching ratios 
compiled in Audi et al. (2003); and (vi) 17 7-ray 
transitions rates for ^^Al (see below). For all nu- 
clear reactions mentioned above, the correspond- 
ing reverse reaction rates were also calculated and 
properly accounted for in the library. Further- 
more, all experimental reaction rates were cor- 
rected for the effects of thermal target excitations 
using the stellar enhancement factors and parti- 
tion functions of Rauscher & Thielemann (2000), 
although it should be noted that these effects are 
relatively small for the sites of nucleosynthesis dis- 
cussed here. For all other reactions for which in- 
sufficient experimental information is available to 
compute reliable experimental rates, the results of 
statistical model calculations (Rauscher & Thiele- 
mann 2000) were adopted. 

The new library, STARLIB, described above ex- 
tends in its present version up to antimony (Sb) 
and is employed for the very first time in the 
present work. A more detailed account of STAR- 
LIB will be published elsewhere (Iliadis et al., in 
preparation). We emphasize that the rate proba- 
bility density functions contained in STARLIB are 
not directly used in the present post-processing 
calculations. However, the tabulated rate factor 
uncertainties will be useful in later sections for the 
discussion of reaction rate uncertainties. 

2.3. Thermal equilibration of ^^Al 

A level scheme of ^^Al is shown in Fig. [TJ The 
5+ ground state, ^^Al^, /3-decays with a half- life of 



f.u. = = \fx'hIgl7hHmi] ^ where Xi denotes a reaction 
rate; this expression can be employed to derive approxi- 
mate rate factor uncertainties, f.u., from published high 
and low reaction rate boundaries. Furthermore, all of the 
replaced rates are obtained from the exact rate tables of 
the original publications, not from any fitting parameters. 



Ti/2 — 7.17 X 10^ y to several excited states (not 
shown in the figure) in the ^^Mg daughter nucleus 
that mainly de-excite via 7-ray transitions to the 
first excited state (at 1809 keV) in 26Mg. The 
subsequent decay to the ^^Mg ground state gives 
rise to the 1809 keV 7-ray line emission observed 
in the Galactic plane (§[!]). An interesting situa- 
tion occurs in ^^Al since its first excited state (228 
keV; 0+), ^^Al™, is an isomer. In other words, the 
significant angular momentum difference between 
ground and isomeric state gives rise to a very small 
7-ray decay constant (with a mean lifetime on the 
order of ~ 10^ y for a M5 decay). Instead, the iso- 
mer /3-decays to the ground state of ^^Mg (without 
emission of a 7-ray) with a half-life of T1/2 = 6.34 
s. 

Although the direct 7-ray transition between 
^^Al^ and 26 jg g^rongly inhibited, they 
may nevertheless be linked (or "communicate") 
with each other via thermal excitations involving 
higher-lying ^^Al levels. Three of these levels are 
shown in Fig. [1] at excitation energies of 417 keV 
(26A1°), 1058 keV (^^Al'') and 2070 keV i^^KY). 
For example, the ground and isomeric state are 
linked via the 7-ray transitions o 417 ^ 228. 
It is obvious that at sufficiently high temperatures 
and long timescales ^^Al^ and ^e^im achieve 
thermal equilibrium, implying that their abun- 
dance ratio is simply determined by the Boltz- 
mann distribution (i.e., by plasma temperature 
and energy difference), and that the internal equi- 
libration mechanism is not important. However, 
it is also clear that, by lowering temperature and 
time scale, the ground and isomeric state will fall 
out of thermal equilibrium at some point. The 
important issue to consider for a reaction net- 
work calculation is if and when ^^Al should be 
treated as a single species (^®A1*, denoting ther- 
mal equilibrium) or as two distinct species (^^Al^ 
and ^^Al™) that /3-decay with their characteristic 
laboratory half lives. This issue was investigated 
by Ward & Fowler (1980), who estabhshed the 
generally accepted procedure: above a limiting 
temperature value of T=0.4 GK, the nuclide ^® Al 
should be regarded as a single species (with all of 
its levels in thermal equilibrium) , while below this 
temperature two distinct species (no equilibrium) 
should be assumed. 

Some of the important 7-ray transition rates in 
^^Al have not been measured yet (indicated by the 
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Fig. 1. — Level scheme of ^^Al. Energies and J'^- 
values are from Endt (1990). The vertical arrows 
represent 7-ray transitions. In ^^Al, the thick ver- 
tical lines denote experimentally measured transi- 
tions (Endt 1990), while the decay rates for the 
thin vertical lines have been estimated using shell 
model calculations (Coc, Porquet & Nowacki 1999; 
Runkle, Champagne & Engel 2001). The possi- 
ble communication of the ^^Al ground state and 
isomeric state (228 keV) via thermal excitations 
involving higher-lying levels (417, 1058 and 2070 
keV) is apparent. The arrows connecting ^® Al and 
^^Mg denote /3-decay transitions, according to the 
shell model calculations of Kajino et al. (1988). 
The /3-decays shown do not necessarily represent 
direct transitions, but rather indicate if the ground 
state or first excited state (at 1809 keV) in ^^Mg 
is predominantly populated by a given transition. 



thin vertical arrows in Fig. [Ij and Ward & Fowler 
(1980) had to employ for these rather crude ap- 
proximations. In more recent work (Coc, Porquet 
& Nowacki 1999; Runkle, Champagne & Engel 
2001) the 7-ray transition rates have been esti- 
mated by using shell model calculations and the 
thermal equilibration of ^^Al was studied in more 
detail using small reaction networks. These works 
confirmed the assumption that above T « 0.45 
GK the ground and isomeric state are in thermal 
equilibrium, while below T 0.15 GK these two 
states /3-decay with their laboratory half lives (i.e., 
neither equilibrium nor any communication be- 
tween levels). In the transitional temperature re- 
gion, T « 0.15 - 0.45 GK, the effective ^^Al decay 
rate was found to deviate substantially from the 
results reported by Ward & Fowler (1980) by up 
to 4 orders-of-magnitude and, furthermore, that 
^^AF and ^e^j™ indirectly linked by 7-ray 
transitions (i.e., they achieve a quasi-equilibrium). 
It was pointed out by Runkle, Champagne & En- 
gel (2001) that these findings rest on simplifying 
assumptions and may not hold if the nuclear reac- 
tions that produce or destroy ^^Al are sufficiently 
fast to disturb the thermal equilibration. See also, 
Gupta & Meyer (2001). 

In the present work we examine carefully the 
equilibration of •^^Al for each of the three nucle- 
osynthesis sites mentioned above. Two separate 
post-processing calculations using recommended 
interaction rates are performed and the resulting 
^^Al yields are compared: one assuming either 
a single or two separate ^^Al species, depending 
on the temperature regime (according to Ward & 
Fowler 1980), and one where the communication 
between ground and isomeric states is explicitly 
taken into account. For the latter case, no arti- 
ficial assumptions about the equilibration of ^^Al 
are made, but additional ^®A1 species (i.e., levels 
at 417, 1058 and 2070 keV; see Fig. [J) need to be 
taken into account in the reaction network. The 
7-ray transition rates linking the various levels in 
^^Al, and the most important /3-decay transition 
rates to ^^Mg, are adopted from Runkle, Cham- 
pagne & Engel (2001). Since these have not been 
published elsewhere, we list the decay constants 
versus temperature in Appendix [X] Furthermore, 
when taking more than two ^^Al species into ac- 
count, the reaction rates for ^^Al^ and '^^AY^ are 
separately needed. This represents an additional 
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complication since, as will be seen below, some of 
these rates are poorly known at present. When 
performing the two separate post-processing cal- 
culations referred to above, it is of utmost impor- 
tance that an internally consistent set of rates for 
the production and destruction of ^^Al', ^^Al^ and 
26^jm employed. Suppose, for example, that the 
rate of a destruction reaction ^^AF(a,b) in the first 
calculation (i.e., assuming a single ^^Al species) is 
adopted from one specific source, and the rates 
of the physically related reactions ^^Al^(a,b) and 
^^Al'"(a,b) in the second calculation (i.e., assum- 
ing five different ^^Al species) are adopted from 
a different source. In such a case, when compar- 
ing the results of the two network calculations, 
one may find significant differences in ^^Al yields. 
However, these may not be caused at all by the 
effects of thermal equilibration of ^^Al but may 
rather reflect spurious results caused by using in- 
consistent reaction rates. In order to elucidate this 
issue, we list in Appendix [B] the rates and sources 
of all reactions considered in our network that pro- 
duce and destroy ^^Al. 

3. Procedure and Results 
3.1. General considerations 

The same reaction network is used in the 
present work for investigating the nucleosynthesis 
of ^^Al in the predicted main locations of massive 
stars: explosive Ne/C burning, convective shell 
C/Ne burning and convective core H burning. 
The network extends from to ''"Ca, including 
175 proton- and neutron-rich nuclides up to the 
respective driplines, that are linked by 1648 nu- 
clear interactions \2.2\i . Initial abundances are 
listed in Tab. [T]and the temperature-density time 
evolutions for each of the sites will be discussed 
in the following subsections. The nucleosynthesis 
will initially be visualized by considering so-called 
"abundance flows", which represent the change 
of abundance per time as a result of an interac- 
tion between two nuclides. Since a forward and 
corresponding reverse reaction occur concurrently, 
what is of main interest is the net abundance flow 
(i.e., the difference between forward and reverse 
flow). Initially, a "standard" network calculation 
is performed, employing recommended reaction 
rates. The final abundances, achieved at the end 
of the standard calculations, are summarized in 



Tab. [H Subsequently, a series of network cal- 
culations is performed, where the rates of many 
reactions are varied individually by generic factors 
of 10, 2, 0.5 and 0.1. Resulting abundance changes 
of ^^Al are then analyzed in detail. Then we fo- 
cus our discussion on the actual rate uncertainties 
in the most relevant temperature region, which 
differs from site to site. It will become appar- 
ent that this temperature region is rather narrow, 
which significantly reduces any (unknown) sys- 
tematic effects that are potentially caused by an 
incorrect temperature dependence of some rates. 
Finally, the issue of thermal equilibration of ^^Al 
is investigated. 

Our strategy regarding which and how many re- 
action rates to vary was as follows. We started by 
considering the net abundance flows (i.e, the dif- 
ference of total abundance flows between a given 
forward and corresponding reverse reaction), in- 
tegrated over the entire duration of a "standard" 
post-processing calculation. These flows will be 
displayed graphically below for each of the three 
burning regimes. All rates of reactions with net 
fiows within 3 ordcrs-of-magnitude of the maxi- 
mum flow were then selected for the variation pro- 
cedure. Obviously, the forward and correspond- 
ing reverse reaction rate need to be multiplied by 
the same variation factor. We added to this list 
all reactions that either destroy or produce ^^Al, 
^^Al and ^^Mg, if these were not taken into ac- 
count already. Furthermore, a number of selected 
other reactions, for example, ^^C(^^C,n)^^Mg and 
2'^Mg(p,7)25Al, were added to the list. We find it 
unlikely that any other reaction not identified by 
the above procedure has a major impact on ^^Al 
nucleosynthesis in massive stars. 

A number of important issues need to be con- 
sidered in detail when performing any reaction 
rate sensitivity study using post-processing cal- 
culations. First, it is assumed that the nuclear 
reaction rates to be varied do not impact the nu- 
clear energy generation. If a given reaction rate 
variation changes both the energy generation and 
the final ^^Al abundance, then this result has no 
obvious meaning. Clearly, in such cases the rate 
should be varied using the full, self-consistent stel- 
lar model. For example, changing the rate of the 
^*'Ne(7,Q!)^®0 reaction, the process that initiates 
Ne burning, influences the ^^Al abundance, al- 
though the effect is small, as will be seen below. 
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Table 1 

Initial and final mass fractions of present post-processing calculations^ 



Nuclide 




X> 










Solar'' 


xNe/C 


C/Ne 


H 


xNe/C 


C/Ne 


H 


iH 






7.0E-01 






1.4E-06 


7.11E-01 


= H 






5.0E-05 








2.76E-05 


^He 






3.1E-05 








3.40E-05 


*He 






2.9E-01 






9.8E-01 


2.74E-01 


12^ 


1.9E-02 


1.5E-01 


3.2E-03 


1.5E-02 


l.OE-01 


2.6E-04 


2.44E-03 


"c 






3.8E-05 






7.3E-05 


2.96E-05 






1.5E-05 






2.7E-06 








8.2E-07 


1.2E-04 


1.2E-03 




2.0E-05 


1.3E-02 


7.90E-04 








4.6E-06 








3.11E-06 




5.0E-01 


7.0E-01 


l.OE-02 


5.8E-01 


6.7E-01 


l.OE-04 


6.55E-03 






2.9E-05 


4.1E-06 




3.7E-06 


4.0E-07 


2.60E-06 


18q 






2.3E-05 








1.48E-05 


""No 


3.4E-01 


l.OE-01 


1.7E-03 


1.8E-01 


1.7E-01 


1.6E-03 


1.16E-03 


"iNc 


2.6E-04 


7.0E-04 


4.3E-06 


3.2E-06 


1.2E-04 




2.92E-06 


""Nc 


1.4E-04 


l.OE-03 


1.4E-04 




l.OE-04 


2.9E-05 


9.37E-05 


""Na 


2.1E-06 


4.2E-06 




1.4E-07 


4.4E-07 








7.5E-03 


9.6E-03 


3.5E-05 


3.5E-04 


8.5E-03 


2.7E-04 


3.84E-05 


^*Mg 


4.6E-02 


4.6E-03 


5.4E-04 


8.7E-02 


1.2E-02 


5.4E-04 


5.61E-04 


"'Mg 


l.OE-02 


5.9E-03 


7.1E-05 


2.2E-03 


6.1E-03 


3.8E-06 


7.40E-05 




8.7E-03 


9.4E-03 


8.1E-05 


2.1E-03 


l.lE-02 


6.0E-05 


8.47E-05 


Al 


o.oE-Uo 


5.8E-07 




1.3E-04 


2.4E-0O 


4.9E-0O 




^''Al 


l.lE-02 


8.6E-03 


6.1E-05 


l.lE-02 


9.5E-03 


l.lE-04 


6.59E-05 




3.7E-02 


l.lE-03 


6.9E-04 


6.7E-02 


2.3E-03 


6.9E-04 


7.49E-04 




1.4E-03 


2.3E-04 


3.6E-05 


1.2E-02 


6.8E-04 


3.6E-05 


3.94E-05 


™Si 


6.7E-04 


1.8E-04 


2.5E-05 


1.2E-02 


3.6E-04 


2.5E-05 


2.69E-05 


31p 


3.9E-04 


8.3E-05 


8.6E-06 


2.5E-03 


l.OE-04 


8.6E-06 


7.53E-06 


32p 


1.5E-06 


1.6E-06 




3.4E-06 


3.2E-06 






33p 


3.3E-07 


7.3E-07 




2.1E-06 


1.4E-06 






32g 


1.4E-02 


2.4E-04 


4.2E-04 


1.6E-02 


2.2E-04 


4.2E-04 


3.92E-04 


33g 


1.9E-04 


2.7E-05 


3.4E-06 


3.8E-04 


5.2E-06 


3.4E-06 


3.20E-06 


34s 


2.7E-03 


l.OE-04 


2.0E-05 


1.9E-03 


l.OE-04 


2.0E-05 


1.85E-05 


35 g 


1.3E-06 


6.9E-07 




3.1E-05 


2.0E-06 






36 


1.9E-06 


9.5E-07 




3.5E-05 


2.2E-06 




7.93E-08 




2.9E-04 


4.9E-06 


2.7E-06 


1.6E-03 


7.1E-06 


2.7E-06 


4.03E-06 


36qj 


7.7E-07 


3.1E-07 




2.1E-05 


2.8E-07 






37 CI 


3.4E-05 


5.1E-05 


9.0E-07 


4.1E-05 


4.9E-05 


9.2E-07 


1.36E-06 


36 Ar 


9.6E-04 


1.7E-05 


8.1E-05 


8.8E-04 


1.3E-05 


8.1E-05 


9.05E-05 


^•^Ar 


1.2E-05 






1.5E-05 








3«Ar 


1.5E-03 


2.7E-05 


1.6E-05 


8.2E-04 


2.9E-05 


1.6E-05 


1.74E-05 


■^»Ar 


2.9E-06 


8.5E-07 




l.lE-05 


2.6E-06 






39j^ 


1.3E-04 


4.7E-06 


3.6E-06 


7.3E-04 


4.2E-06 


3.6E-06 


3.89E-06 



^Only mass fractions of stable or long-lived nuclides in access of X ^ 5 X 10~^ are listed 
here. For ^*^A1 the listed values refer to the ground state. The labels "xNe/C", "C/Ne" and 
"H" refer to explosive Ne/C burning, convective shell C/Ne burning, and convective core H 
burning, respectively. 

'^Initial mass fractions at beginning of post-processing calculations (from Limongi & Chieffi. 
2006). 

*^Final mass fractions at the end of post-processing calculations, obtained from the "stan- 
dard" calculations (i.e., without any reaction rate variations). 

*^Solar system mass fractions, for comparison (from Lodders 2003). 
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Throughout this work, we carefully checked that a 
rate variation did not impact at the same time the 
energy generation. Second, it is very important to 
verify that a given temperature-density-time evo- 
lution is followed precisely in a post-processing cal- 
culation. The time step is numerically adjusted to 
track the abundance evolutions above some lim- 
iting abundance value. For example, spurious 
abundance variations may occur if a time step 
misses the peak temperature even by a few per- 
cent. Therefore, we carefully checked that the 
temperature-density evolution is followed closely 
(within a fraction of a percent). 

Finally, we need to address the issue of stel- 
lar versus laboratory /3-decay rates. The stellar 
plasma affects /3-decays in a number of ways. 
First, at high temperatures thermally excited 
states in the /3-decaying nuclide may undergo tran- 
sitions to levels in the daughter nuclide. Second, 
at high (electron) densities the decay constants 
for electron (or positron) capture will increase. 
Both effects generally cause a change in the total 
/3-decay rate. Many previous stellar model studies 
employed the stellar /3-decay constants calculated 
by Fuller, Fowler & Newman (1982) for the mass 
range of A=21-60, or the more modern results 
of Oda et al. (1994), that are based on shell 
model calculations, for the mass range of A=17- 
39. For the purposes of the present work, the 
highest temperature and density values are en- 
countered in explosive Ne/C burning, with peak 
values of T = 2.3 GK and p = 3.2 x 10^ g/cm^ 
(see below). By inspecting the tables of Oda et 
al. (1994) in the region A<30, we find that for 
the temperatures and densities of interest here 
the stellar and laboratory /3-decay rates are very 
similar in magnitude. The only exceptions are the 
/3-decays of the long-lived species ^^Na and ^^Na. 
However, their destruction via the processes (7,n), 
(7,p) or (p,7) is much faster compared to the /3- 
decays and, therefore, their stellar /3-decay rate is 
unimportant at high values of T and p. In conclu- 
sion, it is sufficient to adopt laboratory /3-decay 
rates throughout this work, except for the /3-decay 
of (thcrmalizcd) ^^Al*. For the calculation of this 
decay we only take the ground and isomeric state 
into account. The decay constants are listed in 
Appendix [X] and agree with the more comprehen- 
sive results of Oda et al. (1994) for temperatures 
and densities below 5 GK and 10^ g/cm^, respec- 



tively. 

3.2. Explosive Ne/C burning 

3.2.1. Standard calculation 

At the beginning of the burning, the most abun- 
dant nucfides are (in order) ^^O, ^°Ne, ^^Mg, ^^Si 
and ^^C (Tab. [T]). The temperature-density- 
time profile for simulating explosive Ne/C burn- 
ing is shown in Fig. [2j It has been extracted 
from a stellar model calculation of a 2QMq star 
(§ 12. ip . Specifically, we select a mass coordinate 
of 2.O4M0, corresponding to the zone where the 
maximum abundance of ^®A1 is produced during 
the explosion. Temperature and density peak at 
T = 2.3 GK and p = 3.2 x 10^ g/cm^, respectively. 
The evolution is followed in a post-processing sim- 
ulation over a total time oi t = 12.8 s. At this 
point the temperature has declined to T = 0.4 GK 
and no additional ^^Al synthesis is occurring. We 
assume at this stage thermal equilibrium for ^^Al, 
i.e., the network contains only a single species, 

26 Al*. 

The net abundance flows, integrated over a 
total running time of t = 12.8 s, for the "stan- 
dard" calculation are displayed in Fig. [3l They 
provide a flrst impression regarding the nucle- 
osynthesis and indicate the degree of "nuclear 
activity". The network consists of all nuclides 
shown as squares. The strongest net abundance 
flows, i.e., those within one, two, and three or- 
ders of magnitude of the maximum flow, are 
displayed by the thickest arrows, arrows of in- 
termediate thickness, and the thinnest arrows, 
respectively. The strongest net flows belong to 
the reactions 20Ne(7,a)i6O and 20Ne(a,7)24Mg 
that drive explosive Ne burning. The released 
a-particles induce a network of secondary re- 
actions, giving rise to a small, but significant, 
abundance of light particles. During the explosive 
burning the mass fractions of protons, a-particles 
and neutrons reach maximum values of 1.0x10"^, 
1.1x10^'°^ and 3.8xl0~^^, respectively. The main 
direct process of ^6^1* synthesis, in terms of 
the net abundance flow, is 25]V[g(p,7)26Al*, with 
2^Mg produced by the '^'^Mg{'D.,^Y^Mg reaction. 
The main neutron sources are 25Mg(a,n)28Si and 
2'5Mg(a,n)29Si. On the other hand, ^^Al* is pre- 
dominantly destroyed via 26Al*(n,p)26Mg and, to 
a lesser degree, 26Al*(n,a)2'^Na. In particular, the 
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Explosive Ne/C burning 




0.01 0.1 1 10 100 



Time (s) 

Fig. 2. — Temperature-density-time evolution for 
explosive Ne/C burning. The profile was ob- 
tained from a hydrodynamic model of a 20Mq star 
(Limongi & Chieffi 2006). In the present post- 
processing calculation, the evolution was only fol- 
lowed until t — 12.8 s (vertical lines), since for 
later times T and p decline to values where no 
additional ^^Al synthesis is occurring. The peak 
temperature and density, near t — 0.6 s, amount 
to T = 2.3 GK and p = 3.2 x 10^ g/cm^, respec- 
tively (see text). 



^^Al*(p,7)^'^Si reaction is entirely negligible un- 
der explosive burning conditions. These general 
features have already been discussed by Limongi 
& Chieffi (2006). The abundance evolutions of 
the species ^^Al* and ^^Al are shown in Fig. 2] 
While the abundance of the latter nuclide is ap- 
proximately constant throughout the calculation, 
the abundance of the former species increases by 
more than an order of magnitude during the ex- 
plosion. The abundance ratio, X (^^ Al*) / X (^"^ Al) , 
increases from an initial value of 3.2 x lO"** to a 
final value of 1.2 x 10"^ (see Tab. [J). 

3.2.2. Reaction rate variations 

Subsequently, the rates of 70 pairs of forward 
and reverse reactions were varied. Those reactions 
whose rate changes have the strongest effect on the 
final ^^Al yield (i.e., at time t = 12.8 s) are listed 
in Tab. [5] All other rate changes, as well as those 
labeled by in the table, produced ^^Al* abun- 
dance changes of less than 20%. The reactions 
are listed in approximate order of importance, as 
measured by their impact on the final ^^Al abun- 
dance. The last two columns display the source 
of the rate and the reported rate uncertainty at a 
temperature near the peak of the explosion (i.e., 
where most of the nucleosynthesis is occurring). 
Disregarding for a moment the actual rate uncer- 
tainties, the six reactions with the strongest im- 
pact on ^^Al* nucleosynthesis are: ^^Al*(n,p)^^Mg, 
25Mg(p,7)26Al*, 25Mg(a,n)28Si, 24Mg(n,7)25Mg, 
20Ne(a,7)24Mg, and 30Si(p,7)3ip. The first and 
second reaction destroys and produces ^^Al*, re- 
spectively, while the third and fourth reaction de- 
stroys and produces, respectively, the ^^Mg seed. 
The fifth reaction produces ^''Mg, from which 
^^Mg is synthesized via neutron capture. The sen- 
sitivity of the final ^^Al* abundance to any of these 
rate changes is thus not surprising. 

The manner by which the sixth reaction im- 
pacts the synthesis of ^^Al* is interesting. In fact, 
the sequence ^°Si(p,7)^^P(p,a)^^Si is the main 
consumer of free protons (together with the pro- 
ton captures on ^®Mg and ^^Al). When the rate 
of the •^*'Si(p,7)^^P reaction is reduced by an or- 
der of magnitude, the number of available pro- 
tons increases near the peak of the explosion and, 
consequently, more ^^Mg nuclei are converted to 
^^Al*. There are 13 more reactions listed in Tab. 
[5] and their mechanisms by which they impact 
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Explosive Ne/C burning 



(a,n) {o:,y) 




Fig. 3. — Net abundance flows, obtained for a post-processing network calculation of explosive Ne/C burning, 
integrated over a total running time oi t — 12.8 s. The T-p profile for this simulation is shown in Fig. [21 
The network consists of all nuclides shown as squares. The strongest net abundance flows, i.e., those within 
one, two, and three orders of magnitude of the maximum flow, are displayed by the thickest arrows, arrows 
of intermediate thickness, and the thinnest arrows, respectively. Thermal equilibrium for ^^Al has been 
assumed (i.e., the network contains only a single species, ^^Al*). 
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Table 2 

Factor changes of final ^^Al* abundance resulting from reaction rate variations for 

EXPLOSIVE Ne/C burning** , ASSUMING THERMAL EQUILIBRIUM FOR ^^Al 



Reaction*^ Rate multiplied by 





10 


2 


0.5 


0.1 


Source'^ 


Uncertainty 


2^Mg(Q;,n)28Si 


0.10 


0.49 


1.8 


4.0 


nacr 


18% 


^^Mg(n,7)2'^Mg 


5.2 


1.6 


0.61 


0.24 


ka02 




0.14 


0.58 


1.6 


3.2 


prGSGIlt 




25Mg(p,7)26Al* 


1.7 


1.4 


0.58 


0.14 


illO 


4% 


30Si(p,7)3^P 


0.51 


0.77 


1.3 


2.0 


illO 


14% 


20Ne(a,7)24Mg 


1.8 


1.4 


0.64 


0.28 


illO 


11% 


27Al(a,p)30Si 


1.5 






0.72 


rath 




29Si(a,n)32s 


0.65 






1.3 


rath 




24Mg(a,7)28Si 


0.62 








iUO 


6% 


24Mg(a,p)27Al 








0.65 


illO 


6% 


27Al(p,7)28Si 


0.60 








illO 


3% 


25Mg(a,p)28Al 
26Al*(n,a)23Na 


0.59 








rath 




0.55 








present 




25Mg(n,7)26Mg 


0.75 








ka02 




28Si(n,7)29Si 


1.4 








ka02 




29Si(p,7)30p 


1.4 








illO 


7% 


32S(n,7)33s 


1.2 








ka02 




26Al*(a,p)29Si 


0.72 








rath 




26Mg(p,7)27Al 


1.2 








illO 


4% 



^The temperature-density-time profile is extracted from a hydrody- 
namic model of a 2OM0 star of initial solar metallicity, see Limongi & 
Chieffi (2006). 

'^In total, the rates of 70 different reactions were varied. Listed are 
only those reactions whose rate changes have the strongest effect on the 
26 Al* yield. All other rate changes, as well as those labeled by 
produced abundance changes of less than 20%. The reactions are listed 
in approximate order of importance. Thermal equilibrium for ^^Al has 
been assumed, i.e., the network contains only a single species, ^^Al*. 

'^Reaction rate references: (nacr) Angulo et al. 1999 (NACRE); (ka02) 
Dillmann et al. (2006) (KADoNiS vO.2); (rath) Rauscher & Thielemann 
(2000); (illO) Iliadis et al. (2010); (present) hybrid rates, see Appendix 
Icelandic:! 

'^Reaction rate uncertainty near a temperature of 2.3 GK, at the peak 
of the explosion; no entry implies that the rate uncertainty is difficult 
to quantify (see text). 
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the final ^^Al* abundance can be easily deduced 
from arguments similar to those given above. The 
only reaction that we found to influence some- 
what the 26 Al* yield but is not listed in Tab. H 
is 2°Ne(7,a)"'^^0. Varying this rate by a factor of 
10 changes the ^^Al* abundance by a factor of ~2. 
However, the estimated rate uncertainty of this re- 
action amounts to only 13% (Iliadis et al. 2010) 
and thus the actual effect is relatively small. 

3.2.3. Reaction rate uncertainties 

Before proceeding, notice the sources of our re- 
action rates, listed in column 6 of Tab. Of 
the 19 reactions listed, the rates of: (i) 8 reac- 
tions are available from the Monte Carlo proce- 
dure (Iliadis et al. 2010; § 12. 2^ : (ii) 6 reactions 
are adopted from the statistical model (Rauscher 
& Thielemann 2000); (iii) 4 reactions are obtained 
from KADoNiS vO.2 (Dillmann et al. 2006); and 
(iv) only one is adopted from NACRE (Angulo et 
al. 1999). Note that none of these rates rely any- 
more on outdated information from Caughlan & 
Fowler (1988). 

We now turn to a discussion of reaction rate 
uncertainties. These are listed for a temperature 
of r — 2.5 GK, near the peak of the explosion, 
in the last column of Tab. [21 when reported in 
the original source. Rate uncertainty estimates 
are of obvious importance. Suppose a rate vari- 
ation of a particular reaction by a factor of 10 
changes the ^^Al* abundance by the same factor. 
Then one may conclude that this particular reac- 
tion rate should be known with rather small un- 
certainty. On the other hand, if a particular rate 
variation barely affects the abundance of ^^Al*, 
one may tolerate a much larger uncertainty. In re- 
ality, however, the issue is much more complicated 
and one is usually confronted with the following 
questions when considering rate uncertainties re- 
ported in the literature. What is the (statistical) 
meaning of a reported rate uncertainty? Is a pre- 
sumed experimental rate at a given temperature 
directly based on data, or is it based on a nor- 
malization of (theoretical) Hauser-Feshbach rates? 
Even if a rate is directly based on data, how large 
is the stellar enhancement factor that must usually 
be obtained from Hauser-Feshbach models? What 
may one estimate for a rate uncertainty if no val- 
ues are reported in the literature? And, even if 
a given rate is directly based on data and if the 



stellar enhancement factor is negligible at a given 
temperature, are there possible systematic errors 
that were not taken into account in the reported 
rate uncertainty? All of these issues play an im- 
portant role and thus reported rate uncertainties 
are frequently difficult to assess. Below we will 
give a few examples to emphasize these points. 

The first reaction listed in Tab. [H '^^Mg{a,nf^Si, 
strongly affects the final ^^Al* abundance. Vary- 
ing the rate by a factor of 10 (2) changes the 
26 Al* yield by a factor of 0.1 (0.5). The rate is 
adopted from the NACRE compilation (Angulo et 
al. 1999), and its reported uncertainty of «18% 
near T = 2.5 GK may on first sight indicate a 
rather reliable rate. However, not enough infor- 
mation is provided in Angulo et al. (1999) to un- 
derstand how exactly this value of uncertainty has 
been obtained. Also, beyond T — 2 GK, i.e., the 
highest temperature for which the rate is directly 
based on data, the rate was extrapolated with 
the aid of (theoretical) Hauser-Feshbach model 
results. Furthermore, even at the lower tempera- 
tures, the rate seems to be based on data from an 
unpublished thesis. Considering these arguments 
together with the importance of the 2^Mg(a,n)2*Si 
reaction, there is no doubt in our minds that this 
particular reaction should be a target of future ex- 
perimental work. Consider, on the other hand the 
fourth reaction listed in Tab. H 25Mg(p,7)26Al*. 
Varying the rate by a factor of 0.1 (0.5) changes 
the 26 Al* yield by a factor of 0.14 (0.58). A rate 
uncertainty of only 4% near T = 2.5 GK has 
been reported by Iliadis et al. (2010). This value 
has been obtained from a Monte Carlo proce- 
dure, implying a statistically meaningful proba- 
bility coverage (68%). The rate near the peak of 
the explosion is directly based on data, i.e., no 
extrapolation using theoretical reaction models is 
needed. Furthermore, the experimental rate is 
normalized to a well-known standard resonance 
strength (for details, see Iliadis 2007). In conclu- 
sion, at present there is less compelling reason for 
remeasuring this reaction at higher energies com- 
pared to the previous case. We emphasize again 
that each reaction must be treated as a special 
case and that a reported rate uncertainty needs to 
be considered carefully. For readers interested in 
the present status of specific reactions, we provide 
brief discussions in Appendix [C] 

The set of rates shown in Tab. [2] that are based 
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on Hauser-Feshbach theory (labelled by "rath") 
represent a special case. It has been claimed by 
Rauscher & Thielemann (2000) that "...the accu- 
racy of the rates is estimated to be within a factor 
of 1.5-2...". Obviously, if too few resonances con- 
tribute to the rate at a given temperature, the 
statistical model will provide a poor description. 
For this reason, Rauscher & Thielemann (2000) 
provide a minimum temperature estimate below 
which the Hauser-Feshbach rates become inaccu- 
rate. This minimum temperature value is calcu- 
lated from a parameterization of nuclear level den- 
sities, assuming that at least 10 levels (Rauscher, 
Thielemann & Kratz 1997) are located in the as- 
trophysically important energy window (e.g., the 
Gamow peak for charged-particle reactions). Note 
that for all of the reactions labeled "rath" in 
Tab. [2] the peak temperature of the explosion 
(T = 2.3 GK) far exceeds the minimum tempera- 
ture required for the applicability of the Hauser- 
Feshbach model according to Rauscher & Thiele- 
mann (2000). Unfortunately, the above claims 
are not supported when comparing the Hauser- 
Feshbach rates with results that are directly based 
on experiment. The issue was discussed in Iliadis 
et al. (2001), who found that for several reactions 
involving A=20-40 mass targets "... the deviation 
between theoretical and experimental rates far ex- 
ceeds the usually quoted factor of 2 reliability of 
statistical model results ..." . Clearly, more work is 
required to resolve this controversy. At this point 
it may be argued that all of the reactions labeled 
by "rath" in Tab. [2] should be targets for future ex- 
perimental work, including the important destruc- 
tion reactions 26Al*(n,p)26Mg and ^^Al\n,a)^^'Nii. 
(labeled "present"; see Appendix lC.6l and lC.5p . 

Neutron capture rates represent another spe- 
cial case. We adopted for these the results pre- 
sented in the KADoNiS vO.2 evaluation (Dillmann 
et al. 2006; these rates are labelled by "ka02" in 
Tab. [2|). The most important neutron capture 
reaction for the purposes of the present work is 
^'*Mg(n,7)^^Mg, as is apparent from the table. In 
order to obtain a better sense for the uncertainties, 
we will briefly discuss how the KADoNiS evaluated 
rates have been obtained and what information 
is actually incorporated in reaction rate libraries. 
The arguments below apply equally to the other 
reactions listed in Tab. [H i.e., 25Mg(n,7)2'5Mg, 
28Si(n,7)29Si and ^^S{n,-f)^^S. The KADoNiS eval- 



uation tabulates recommended rates for the range 
of A;r = 5 - 100 keV (corresponding to T = 0.06 - 
1.2 GK). For the neutron captures on '^'^'^^Mg, ^^Si 
and "^^S the rates are obtained from experimental 
data on resonance properties (with some theoret- 
ical corrections for direct neutron capture contri- 
butions, if applicable) over the entire tabulated 
temperature range. According to the KADoNiS 
evaluation, the "relative uncertainties [of the rates] 
are similar to those quoted for the 30 keV data" 
(12% for neutron capture on ^^Mg). The tabu- 
lated rates include the stellar enhancement factor 
(§ 12. 2p . although these are predicted to be close 
to unity for the (n,7) reactions mentioned above. 
Note that for explosive Ne/C burning the rates 
are needed at temperatures (T w 2.3 GK) that 
have not been covered by experiments. Thus it 
is not obvious how to extrapolate the rates from 
lower temperatures, where they are based on ex- 
perimental data, to much higher temperatures. 
Furthermore, it must be pointed out that the ex- 
perimental KADoNiS rates are not directly used 
in reaction rate libraries (including ours). What 
is usually incorporated for neutron captures are 
Hauser-Feshbach rates, which are normalized to 
the experimental rates at a single temperature 
{kT = 30 keV or T = 0.35 GK). Since the level 
density for targets in the mass A < AO range at 
kT = 30 keV may be too small for the applica- 
tion of statistical models, an additional system- 
atic uncertainty is introduced when extrapolating 
such normalized rates to higher temperatures. For 
example, in the case of ^^Mg(n,7)^^Mg the (exper- 
imental) KADoNiS rate at the upper temperature 
cutoff {T = 1.2 GK) deviates from the normalized 
Hauser-Feshbach rate already by «40%. Consider- 
ing the above arguments, we estimate a rate uncer- 
tainty of a factor of 2 — 3 for the neutron captures 
on ^^'^^Mg, ^*Si and '^^S near the peak of explosive 
Ne/C burning. Recall from Tab. [2] that varying 
the ^'*Mg(n,7)^^Mg rate by a factor of 2 increases 
the ^^Al* yield by a factor of 1.6. Clearly, a more 
reliable experimental rate for ^*Mg(n,7)^^Mg at 
higher temperatures is urgently needed. 

3.2.4- Thermal equilibration 

We will now consider the issue of thermal equi- 
libration. Recall that we assumed so far a sin- 
gle species of ^^Al, implying thermal equilibrium 
(^^Al'). We will now relax this assumption and 
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follow the equilibration numerically in the network 
calculation. To this end, we introduce five differ- 
ent species of •^^Al, as explained in § 12.31 The 
required 7- and /3-decay transitions between and 
from these levels are discussed in detail in Ap- 
pendix |X1 The price we pay is that additional re- 
action rates, involving ^^Al^ and '^^AV" separately, 
have to be incorporated into the network (see Ap- 
pendix [B|). As will be seen, some of these rates are 
highly uncertain. 

As a first step, we performed a standard 
post-processing network calculation (with rec- 
ommended rates) using the same temperature- 
density-time evolution as before (Fig. [2]). The 
final ^^Al^ abundance, at t = 12.8 s, is found to 
be identical to our earlier result obtained assum- 
ing thermal equilibrium (Tab. [IJ. Thus the latter 
assumption seems to be justified. An impression 
can be gained from Fig. [SJ showing the abundance 
evolution of ^^Al levels. The top part displays the 
abundances for individual ^^Al species and it is 
apparent that at any given time the ^^AF abun- 
dance dominates over those of the other species. 
The bottom part displays the fraction of the to- 
tal ^^Al abundance that resides in the isomeric 
state. This curve is directly obtained from the 
network calculation, but is indistinguishable from 
the one calculated assuming a Boltzmann distri- 
bution (i.e., thermal equilibrium). 

Next, the rates of 20 different reactions and 
transitions, together with their inverse rates, were 
varied individually by factors of 100, 10, 2, 0.5, 
0.1 and 0.01. This list contained all nuclear reac- 
tions that produced and destroyed ^^Al" or ^^Al™. 
It also included those f3- and 7-ray decay rates of 
^^AP levels that were estimated using the shell 
model, as described in § 12.31 and shown in Fig. 
[TJ Experimentally obtained /3- and 7-ray decay 
rates have not been varied since their uncertain- 
ties are very small. The final ^^Al^ abundance, 
after each rate variation, was then compared to 
the standard calculation. The results are listed in 
Tab. [3] It can be seen that of the 20 interactions 
only eight, all of them reactions, influence the fi- 
nal ^^AV yield. In other words, even a variation 
by a factor of 100 of the shell-model based (3- or 
7-ray decay rates seems to have no effect on the 
^^Al^ abundance. The eight reactions displayed 
in Tab. [3] are listed in approximate order of im- 
portance. It must be noted that the impact of 



these rates, with one exception, on the final abun- 
dance of ^^AW is moderate. For example, consider 
the second and third reaction, ^^Mg(p,7)^^AF and 
^^Mg(p,7)^^Al'". Even a rather small rate varia- 
tion (by a factor of 2) infiuences the final ^^Al 
yield (by «30%). However, these rates are based 
on experimental information (Iliadis et al. 2010) 
and their Monte Carlo uncertainties are predicted 
to amount to only 4-5% near the peak of the ex- 
plosion. 

The one exception is the ^^AP(n,p)^^Mg reac- 
tion rate. Increasing this rate by a factor of 100 
changes the ^^Al^ abundance by a factor of 0.017. 
We adopted for this reaction the same rates as for 
^^Al*(n,p)^^Mg (see Appendix |B|. Our numerical 
results indicate that even a factor of 100 varia- 
tion in the ^^AF (n,p)^^Mg rate does not change 
the thermal equilibrium abundance ratio of '^^A\™ 
and ^^AW (Fig. [S])- In fact, comparison of the first 
entry of Tab. [3] with the third entry of Tab. [5] im- 
mediately reveals that the ^^Al(n,p)^^Mg reaction 
impacts the final ^^Al abundance by the same fac- 
tor changes, no matter if a single (thermalized) or 
five species of ^^Al are used in the simulation. In 
conclusion, ^^Al is in thermal equilibriurrQ during 
explosive Ne/C burning and, consequently, there 
is no need to introduce the extra complication of 
five ^^Al levels, and their mutual interactions, into 
the reaction network. 

3.3. Convective shell C/Ne burning 

3.3.1. Standard calculation 

Preliminary studies of the impact of nuclear 
uncertainties on pre-explosive ^^Al yields can be 
found in Baldovin, Pignatari & Gallino (2006). 
These authors performed post-processing studies 
using a schematic one-zone model consisting of 
two phases: a constant temperature of T w 1.1 
GK until the ^^C mass fraction decreases from an 
inital value of 0.18 to 0.10 for phase 1, and a con- 
stant temperature of T w 1.3 GK until the ^■^C 
mass fraction reaches a value of 0.050 for phase 2. 



*The reader may suspect circular reasoning in our argu- 
ments, in the sense that we used the same rates for the 
ground and isomeric states as for the thermahzed ^®A1 tar- 
get, and then conclude that ^®A1 is in thermal equilibrium. 
However, this is not the case since our assumption for the 
nominal rates serves as a starting point only and we fully 
explore individual rate changes by factors up to 100. For 
more information, see App. [B] 
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Table 3 

Factor changes of final ^^Al^ abundance resulting from reaction rate variations for 

EXPLOSIVE Ne/C burning'^ , ASSUMING FIVE SPECIES OF ^^Al 



Reaction*^ Rate multiplied by 





100 


10 


2 


0.5 0.1 


0.01 


Source'^ 


Uncertainty 


26AP(n,p)26Mg 


0.017 


0.14 


0.57 


1.6 2.9 


3.8 


present 




25Mg(p,7)26A19 


1.7 


1.6 


1.3 


0.71 0.45 


0.38 


iUO 


4% 


25Mg(p,7)26Al" 


1.6 


1.6 




• • • 0.79 


0.79 


iUO 


5% 


26A19(a,p)29Si 


0.21 


0.71 








rath 




26A19(n,a)23Na 


0.21 


0.54 








present 




26Al"(n,p)26Mg 


0.36 










present 




26AF(p,7)27si 


0.52 










illO 


7% 


26Ar(n,a)23Na 


0.79 










present 





"^The temperature-density-time profile is extracted from a hydrodynamic model of a 
20Mo star of initial solar metallicity, see Limongi & ChiefR (2006). 

''In total, the rates of 20 different reactions producing or destroying ^^Al^ and ^e^jm 
were varied. Listed are only those reactions whose rate changes have the strongest effect 
on the ^^Al^ yield. All other rate changes, as well as those labeled by produced 
abundance changes of less than 20%. The reactions are listed in approximate order 
of importance. No thermal equilibrium for ^^Al has been explicitly assumed, i.e., the 
network contains five different species {^^kW, 26A1™, 26^^^ '^^ ^ '^^kV) and takes the 
interactions between them into account. 

■^Reaction rate references: (illO) Iliadis et al. 2010; (rath) Rauscher & Thielemann 
2000; (present) hybrid rate, see Appendix IC.5I and IC.61 In the latter three cases, we 
assumed that the rate involving ^^Al^ or ^e^jm jg ^j^g game as the rate for ^^Al* (see 
comments in Appendix IbI) . 

''Reaction rate uncertainty near a temperature of 2.3 GK, at the peak of the explosion; 
no entry implies that the rate uncertainty is difficult to quantify. 



17 



In the present work we proceeded as follows. Ini- 
tially, we extracted the temperature-density-time 
evolution of the deepest and hottest zone of the 
convective C/Ne burning shell from a stellar evo- 
lution model of a 60 Mq star with initial solar 
mctallicity (Limongi & ChiefH 2006). This pro- 
file extended in time from the formation of the 
shell until a time of t = 3.15 x 10^ s. Using this 
T — p profile directly in a post-processing study 
would greatly distort the nucleosynthesis predic- 
tion, since the effects of convection are not taken 
properly into account. On the one hand, convec- 
tion constantly carries fresh fuel (here ^^C) into 
the burning region, while, on the other hand, it 
transports fragile nuclei from the burning region to 
cooler layers where they survive for a longer period 
of time. Therefore, in a stellar evolution calcula- 
tion, convection has the effect of lengthening con- 
siderably the duration of nuclear burning. If we 
would use this profile directly in a post-processing 
simulation, then the initial ^^C fuel, for example, 
would be destroyed much faster than in the actual 
stellar evolution calculation. After some trial at- 
tempts, we found that compressing the time axis 
of the original T — p profile by a factor of 60 gives 
results that are consistent with the stellar evolu- 
tion calculations. Clearly, this large scaling fac- 
tor reflects the strong effects of convection during 
C/Ne shell burning. The results are shown in Fig. 
m displaying the temperature and density depen- 
dence on the ^^C mass fraction as solid and dashed 
lines, respectively. For comparison, the circles in- 
dicate the corresponding values from the stellar 
evolution calculations. The good agreement is en- 
couraging and thus we used the scaled T — p profile 
for our post-processing study. At the beginning of 
the burning, when Xi{^'^C) = 0.15, temperature 
and density start at values of T — 1.13 GK and 
p = 6.3 X 10'* g/cm"^, respectively. The profile ex- 
tends over a time period of < = 5.24 x 10^ s, when 
Xfi^'^C) = 0.10, and ends with values of T = 1.44 
GK and p = 1.1 x 10^ g/cm^. 

At the beginning of the burning, the most abun- 
dant nuclides are (in order) -'^^O, ^^C and ^"^Ne 
(Tab. [T|). We assume at this stage thermal equi- 
librium for ^^Al, i.e., the network contains only a 
single species, ^^Al*. The net abundance flows, in- 
tegrated over a total running time oil = 5.24 x 10** 
s, for the standard calculation are displayed in 
Fig. [71 The strongest net flows belong to the -pri- 



mary carbon burning reactions ^^C('^^C,a)^''Ne 
and ^^C(^^C,p)^'^Na, and to the secondary reac- 
tions i'^0(a,7)^°Ne and 23Na(p,a)2°Ne that are 
initiated by the released light particles from the 
primary reactions. Near the end of convective 
shell C/Ne burning the mass fractions of pro- 
tons, a-particles and neutrons reach maximum 
values of 9.5x10"", 5.0x10-^ and 2.9x10-^6, 
respectively. The main direct process of ■^^Al* 
synthesis, in terms of the net abundance flow, 
is 25Mg(p,7)26Al*, with ^^Mg produced by the 
22Ne(a,n)25Mg and 2**Mg(n,7)25Mg reactions. 
The main neutron source is ^^Ne(a,n)^^Mg. On 
the other hand, ^^Al* is mainly destroyed via the 
/?-decay ^^Al* -j>26Mg (see column 2 of Tab. (5]). In 
particular, the neutron abundance is too low in the 
standard calculation for the destruction reactions 
26Al*(n,p)26Mg and 26Al*(n,a)2^Na to compete 
successfully with the /3-decay of ^^Al*. The abun- 
dance evolutions of the species ^^Al* and ^^Al are 
shown in Fig. [51 While the abundance of the latter 
nuclide is approximately constant throughout the 
calculation, the abundance of the former species 
increases by more than an order of magnitude over 
the course of the burning. The abundance ratio, 
X(^^ AV-)IX(^'^ AV), increases from an initial value 
of 6.7 X 10^^ to a final value of 2.5 x 10~^ (see 
Tab. [1]). 

3.3.2. Reaction rate variations 

Subsequently, the rates of 66 pairs of forward 
and reverse reactions were varied. Those reactions 
whose rate changes have the strongest effect on the 
final ^^Al yield (i.e., at the end of the calculation, 
when X/(i2C) = 0.10) are hsted in Tab. H All 
other rate changes, as well as those labeled by 
in the table, produced ^^Al* abundance changes of 
less than 20%. The reactions are listed in approx- 
imate order of importance, as measured by their 
impact on the final ^^Al abundance. The last two 
columns display the source of the rate and the re- 
ported rate uncertainty at a temperature of «1.4 
GK near the end of the calculation. Disregarding 
at first the actual rate uncertainties, the four re- 
actions with the strongest impact on ^^Al* nucle- 
osynthesis are: '^^Mg{p,"ff^Al\ 26Al*(n,p)26Mg, 
23Na(p,a)20Ne and '^^Na{a,p)'^^Mg. The first re- 
action produces ^^Al*, while multiplying the rate 
of the second reaction by a factor of 10 would make 
it the dominant ^^Al* destruction process, at the 
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Convective shell C/Ne burning 



(a,n) (c<,7) 




n 

1 



Fig. 7. — Net abundance flows, obtained for a post-processing network calculation of convective shell C/Ne 
burning, integrated over a total running time of t = 5.2 x 10'' s, when the ^^C mass fraction has decreased 
to 0.097. The T-p profile for this simulation is shown in Fig. [B) The network consists of all nuclides shown 
as squares. The strongest net abundance flows, i.e., those within one, two, and three orders of magnitude of 
the maximum flow, are displayed by the thickest arrows, arrows of intermediate thickness, and the thinnest 
arrows, respectively. Thermal equilibrium for ^^Al has been assumed (i.e., the network contains only a single 
species, ^^Al*). 
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cost of the /?-decay of ^^Al*. The third reaction is 
the main consumer of free protons. When the rate 
of the ^'^Na(p,a)^°Ne reaction is increased, the 
number of available protons decreases and, con- 
sequently, fewer ^^Mg nuclei can be converted to 
^^Al*. The fourth reaction represents the second 
most important proton-generating process (after 
the primary ^^C(^^C,p)^^Na reaction). When the 
^^Na(a,p)^^Mg reaction rate is increased, more 
protons are available for producing ^^Al* from 
25Mg. 

Other important rate variations that im- 
pact the final ^^Al* abundance arise from the 
reactions ^^Mg(a,n)2*Si, 2''Mg(n,7)2^Mg and 
^^Al*(n,Q!)^'^Na, which have already been dis- 
cussed in § 13.2.21 In addition, the reactions 
i6O(a,7)20Ne, i2c(i2c,n)23Mg and 26Mg(a,n)29Si 
play an important role. Decreasing the rate of the 
first reaction will consume fewer a-particles, thus 
increasing the production of neutrons (for con- 
verting ^"^Mg to ^^Mg seed) via ^^Ne(a,n)^'^Mg. 
Increasing the rate of the second reaction pro- 
duces an additional burst of neutrons towards the 
very end of the burning, thereby destroying more 
^^Al* nuclei via the (n,p) and (n,Q:) reactions. For 
specific comments on the ^^C(^^C,n)^3Mg reac- 
tion, see Appendix IC.ll Similar arguments apply 
to the ^^Mg(Q;,n)^^Si reaction. There are four 
more reactions hsted in Tab. IH ^''Al(n,7)^*Al, 
25Mg(n,7)26Mg, 26Mg(p,7)27Al and '^'^ A\{p,a)^'^Mg, 
which impact the final ^^Al* abundance. The only 
reactions that we found to infiuence the Al' yield 
but are not listed in Tab. Hare ^^Ci^^C,^)^^^^ 
and ^^C(^^C,a)^'^Ne. Varying these rates by a 
factor of 10 changes the ^^Al* abundance by a fac- 
tor of Ri2. However, these reactions drive carbon 
burning and thus strongly infiuence the nuclear 
energy generation. Therefore, varying this rate in 
a post-processing study is not very meaningful. 
Nevertheless, the effect appears to be relatively 
small. 

3.3.3. Reaction rate uncertainties 

Of the 14 reactions listed in Tab. |4j the rates 
of: (i) 5 reactions are available from the Monte 
Carlo procedure (Iliadis et al. 2010; §[2J]); (n) 3 
reactions are adopted from the statistical model 
(Rauscher & Thielemann 2000); (iii) 3 reactions 
are obtained from KADoNiS vO.2 (Dillmann et 
al. 2006); and (iv) 2 reactions are adopted from 



NACRE (Angulo et al. 1999). Only the rates of 
the ^^C(^^C,n)^'^Mg reaction are partially based 
(i.e., the total ^^C-|-^^C rate) on the information 
provided in Caughlan & Fowler (1988), see Ap- 
pendix IC.ll 

Reaction rate uncertainties are listed for a tem- 
perature of T = 1.4 GK, near the end of the 
burning, in the last column of Tab. SI when re- 
ported in the original source. The uncertainties for 
23Na(p,a)20Ne, 25Mg(p,7)26Al*, ^^0{a,-if'^'tie, 
2'^Mg(p,7)2'i'Al and 2'^Al(p,a)2''Mg are relatively 
small and, therefore, the rate estimates for these 
reactions seem suflaciently reliable at present. The 
^^Mg(a,n)2^Si and ^^Mg(a,n)2*Si reactions are 
listed with rather large rate uncertainties (29% 
and 59%, respectively, according to Angulo et 
al. 1999) and, therefore, should be addressed in 
future work (see also Appendix IC.3p . No rate 
uncertainties are given for any of the other re- 
actions listed in the table. These rates are de- 
rived, for example, from Hauser-Feshbach the- 
ory (Rauscher & Thielemann 2000) or from the 
KADoNiS vO.2 evaluation (Dillmann et al. 2006) 
and uncertainties are difficult to quantify, as has 
already been discussed in ij |3.2.3l Clearly, more re- 
liable experimental rates for these reactions are ur- 
gently needed. Specific comments on the reactions 
26Al*(n,p)26Mg, 26Al*(n,a)23Na, ^'^^&{a,pf^Mg 
and ^2C(^^C,n)23Mg can be found in AppendixO 

3.3.4. Thermal equilibration 

So far we assumed a single species of ^^Al, im- 
plying thermal equilibrium (^^Al*). We will now 
follow the equilibration numerically in the network 
calculation. Five different species of ^^Al are in- 
corporated into the network, as explained in § 12.31 
The required 7- and /3-decay transitions between 
and from these levels are discussed in detail in Ap- 
pendixlAl The additional reaction rates, involving 
^^AV and 26^1™ separately, are discussed in Ap- 
pendix |B] 

As a first step, a standard post-processing net- 
work calculation (with recommended rates) is per- 
formed using the same temperature-density-time 
evolution as before (Fig. [5]) . The final ^^AV abun- 
dance, when the ^^C mass fraction has fallen to a 
value of Xf{^'^C) = 0.10, is found to be identical 
to our earlier result obtained assuming thermal 
equilibrium (Tab. [1]). Thus the latter assumption 
seems to be justified. An impression can be gained 
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Table 4 

Factor changes of final ^^Al* abundance resulting from reaction rate variations for 

CONVECTIVE shell C/NE BURNING** , ASSUMING THERMAL EQUILIBRIUM FOR ^^Al 



Reaction*^ Rate multiplied by 





10 


2 


0.5 


0.1 


Source'^ 


Uncertainty 


23Na(p,a)2"Ne 
26Al*(n,p)26Mg 


0.15 


0.61 


1.6 


4.2 


illO 


6% 


0.16 


0.65 


1.4 


1.9 


present 




25Mg(p,7)26Al* 


6.2 


2.0 


0.46 


0.10 


iUO 


5% 


23Na(a,p)26Mg 


3.0 


1.3 




0.71 


rath 




26Mg(a,n)29Si 


0.40 


0.83 




1.3 


nacr 


29% 


2'^Mg(n,7)25Mg 
26Al*(n,a)23Na 


2.1 


1.3 




0.70 


ka02 




0.54 


0.79 






present 




i6O(a,7)20Ne 




0.83 


1.3 


1.7 


iUO 


14% 


25Mg(a,n)28Si 


0.42 








nacr 


59% 


i2C(i2C,n)23Mg 


0.46 








da77 




27AI(n,7)28Al 


1.7 








ka02 




2SMg(n,7)26Mg 


1.3 








ka02 




26Mg(p,7)27Al 


0.71 








illO 


5% 


27Al(p,a)24Mg 


0.79 








illO 


7% 



^The temperature-density-time profile is extracted from a stellar evolu- 
tion calculation of a 60AfQ star with initial solar metallicity, see Limongi 
& Chiefs (2006). 

''In total, the rates of 66 different reactions were varied. Listed are 
only those reactions whose rate changes have the strongest effect on the 
26 Al* yield. All other rate changes, as well as those labeled by 
produced abundance changes of less than 20%. The reactions are listed 
in approximate order of importance. Thermal equilibrium for ^^Al has 
been assumed, i.e., the network contains only a single species, ^^Al*. 

'^Reaction rate references: (nacr) Angulo et al. 1999 (NACRE); (ka02) 
Dillmann et al. (2006) (KADoNiS vO.2); (rath) Rauscher & Thielemann 
(2000); (illO) Iliadis et al. (2010); (present) hybrid rate, see Appendix 
[aHandESl (da77) total ^^C+^^C rate from Caughlan & Fowler (1988), 
with neutron branching ratio adopted from Dayras et al. (1977), see 
Appendix IC.ll 

Reaction rate uncertainty near a temperature of 1.4 GK, at the end 
of the calculation; no entry implies that the rate uncertainty is difficult 
to quantify (see text). 
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from Fig. [9l showing the abundance evolution of 
^^Al levels. The top part displays the abundances 
for individual ^^Al species and it is apparent that 
at any given time the ^^AF abundance dominates 
over those of the other species. The bottom part 
displays the fraction of the total ^^Al abundance 
that resides in the isomeric state. This curve is di- 
rectly obtained from the network calculation, but 
is indistinguishable from the one calculated assum- 
ing a Boltzmann distribution (i.e., thermal equi- 
librium) . 

Subsequently, the rates of 20 different interac- 
tions, together with their inverse processes, were 
varied individually by factors of 100, 10, 2, 0.5, 
0.1 and 0.01. This list contained all nuclear reac- 
tions that produced and destroyed ^^Al^ or ^^Al™. 
It also included those /3- and 7-ray decay rates of 
^^AP levels that were estimated using the shell 
model, as described in § 12.31 and shown in Fig. [T] 
Experimentally obtained /?- and 7-ray decay rates 
have not been varied since their uncertainties are 
very small. The final ^^Al^ abundance, after each 
rate variation, was then compared to the standard 
calculation. The results are listed in Tab. El It 
can be seen that of the 20 interactions only five, all 
of them reactions, infiuence the final ^^AF yield. 
In other words, even a variation by a factor of 100 
of the shell-model based /?- and 7-ray decay rates 
seems to have no effect on the ^^Al^ abundance. 

The five reactions displayed in Tab. [S]are listed 
in approximate order of importance. For example, 
increasing the ^^Mg(p,7)^^A15 rate by a factor of 
2 will enhance the final ^^Al yield by w50%. How- 
ever, this rate is based on experimental informa- 
tion (Iliadis et al. 2010) and their Monte Carlo un- 
certainty is predicted to amount to only 5% near 
T — 1.4 GK. As was the case for explosive Ne/C 
burning (§ 13.2.4^ . we found that even a factor of 
100 variation in the rates of these five reactions has 
no impact on the thermal equilibrium abundance 
ratio of ^SAl™ and ^^Alf (Fig. E]). Comparison of 
the factor changes listed in Tab. [5] with those of 
Tab. |4]reveals that these five reactions impact the 
final ^^Al abundance by similar amounts, no mat- 
ter if a single (thermalized) or five species of ^^Al 
are used in the simulation. In conclusion, ^^Al 
is in thermal equilibrium during convective shell 
C/Ne burning and, therefore, there is no need to 
introduce the extra complication of five ^^Al lev- 
els, and their mutual interactions, into the reac- 



tion network. 

3.4. Convective core H burning 

3.4-1- Standard calculation 

We have performed post-processing calcula- 
tions using the temperature-density profile for 
convective core H-burning in an 80 M© star of so- 
lar initial composition (Limongi & Chieffi 2006). 
Following the procedure adopted for convective 
shell C/Ne burning, we have artificially shortened 
the burning time in our calculation so that the 
time evolution of the hydrogen fuel would closely 
follow that for the stellar evolution calculation. 
For this profile, compressing the time axis by a 
factor of 17 gave consistent results, as shown in 
Fig. \M Burning starts at X.^-H) = 0.70, T = 
0.044 GK and p = 2.03 g/cm^. The modified pro- 
file extends for t — 5.85 x 10^^ s, at which time 
XfC^H) = 1.4 X 10-6, T = 0.088 GK and p = 
17.9 g/cm"^. Note, that in this case T and p refer 
to the values at the center of the star. 

Our standard calculation assumes that ^^AF 
and 26j\im a,re distinct species. The net abundance 
flows are shown in Fig. [TT] and not surprisingly, 
the strongest are within the CNO cycles. Both 
the ground state of ^®A1 and the isomeric level 
are produced via the ^^Mg(p,7)^^Al reaction. For 
most of the burning period, ^^Al is produced from 
the reservoir of initial ^^Mg and it is not until the 
very late stages of burning that ^^Mg is replen- 
ished through ^^Mg(p,7)^^Al and the subsequent 
/3-decay of ^^Al. The primary destruction route 
for the isomer is /3-decay to ^®Mg, whereas the 
ground state is destroyed via ^^Als(p,7)^'''Si. The 
abundance evolution of ^^Al^ and ^^Al is shown 
in Fig. [T2I The abundance of ^^Al is essentially 
constant at X(^'''A1) = 6.1 x 10~^ until late times 
while ^^AF grows to a maximum of X(^^AF) = 
5.2 X 10-^ before dropping to Xi^'^AW) = 4.9 x 
10^^ at the end of burning (see Tab. 1). The ra- 
tio ^^AF/^^Al shows a similar behavior, reaching 
a maximum of 0.73 with a final value of 0.45. 

3-4-2- Reaction rate variations, thermal equili- 
bration and uncertainties 

The rates of 26 pairs of forward and reverse 
reactions were varied and those reactions whose 
rate changes have the strongest effect on the final 
abundance of ^^Al^ (i.e., at the end of the calcu- 
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Table 5 

Factor changes of final ^^Al^ abundance resulting from reaction rate variations for 

CONVECTIVE SHELL C/Ne BURNING^ , ASSUMING FIVE SPECIES OF ^^Al 



Reaction*^ Rate multiplied by 





100 


10 


2 


0.5 0.1 


0.01 


Source'^ 


Uncertainty 


26AF(n,p)26Mg 


0.017 


0.16 


0.63 


1.3 1.9 


2.0 


present 




25Mg(p,7)26A19 


2.9 


5.4 


1.5 


0.63 0.35 


0.29 


iUO 


5% 


25Mg(p,7)26Al" 


6.7 


3.0 




• • • 0.75 


0.71 


iUO 


6% 


26A19(n,a)23Na 


0.12 


0.54 








present 




26Al™(n,p)26Mg 


0.58 










present 





^The temperature-density-time profile is extracted from a stellar evolution calculation 
of a 6OM0 star with initial solar metallicity, see Limongi & Chieffi (2006). 

''In total, the rates of 20 different reactions producing or destroying ^^Al^ and ^e^l™ 
were varied. Listed are only those reactions whose rate changes have the strongest effect 
on the ^^Al^ yield. All other rate changes, as well as those labeled by produced 
abundance changes of less than 20%. The reactions are listed in approximate order 
of importance. No thermal equilibrium for ^^Al has been explicitly assumed, i.e., the 
network contains five different species {^^kW, 26A1™, 26 AP, 26^16^ 26^^) ^nd takes the 
interactions between them into account. 

■^Reaction rate references: (illO) Iliadis et al. 2010; (present) hybrid rate, see Appendix 
IC.5l and lC.61 In the latter two cases, we assumed that the rate involving "^^AW or 26^1™ 
is the same as the rate for 26 Al* (see comments in Appendix [B|) . 

Reaction rate uncertainty near a temperature of 1.4 GK, at the end of the calculation; 
no entry implies that the rate uncertainty is difficult to quantify. 
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Core H burning 



Ca,n) [a.y) 




Fig. 11. — Net abundance flows, obtained for a post-processing network calculation of convective core H 
burning, integrated over a total running time of t = 5.9 x 10^^ s, when the mass fraction has decreased to 
1.3 X 10^^. The T-p profile for this simulation is shown in Fig. [TUl The network consists of all nuclides shown 
as squares. The strongest net abundance flows, i.e., those within two, four, and six orders of magnitude of 
the maximum flow, are displayed by the thickest arrows, arrows of intermediate thickness, and the thinnest 
arrows, respectively. 
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lation, when Xh = 1-3 x 10~^) are listed in Tab. 
[6l All other rate changes, as well as those labeled 
by in the table, changed the ^^Al^ abundance 
by less than 20%. The reactions are listed in ap- 
proximate order of importance, as measured by 
their impact on the final ^^Al^ abundance. The 
last two columns display the source of the rate 
and the reported rate uncertainty at a temper- 
ature of w 0.09 GK, near the end of the calcu- 
lation. For these calculations, ^^Al^ and 26^1™ 
are considered to be separate species. It is not 
surprising that the 2^Mg(p,7)2^Als reaction has 
the largest impact on the final ^^Al^ abundance, 
but it is interesting that variations in the rate of 
25Mg(p,7)26Al'" also affect the ^'^AW abundance. 
This is because these two reactions are the only 
significant destruction mechanisms for ^^Mg and 
thus 2'^Mg(p,7)2^Al™ lowers the amount of ^^Mg 
available to be converted to ^^AF. The third re- 
action, 26A19(p,7)27s i, is the major destruction 
route for ^'^AW. Finally, the ^'^0{p,"f)^'^F reaction 
affects ^^Al production by reducing the abundance 
of free protons. At no point during the calculation 
did the temperature reach a point where ^^Al" and 
26 j^Yin f^Q-^i^ communicate through thermal excita- 
tions. This was verified by including the three me- 
diating levels discussed in § 2.3 in a second series 
of network calculations and no change was seen in 
the final abundance of ^^Al^. 

The rates for all of the reactions listed in Tab. 
[S]were obtained using the Monte Carlo procedure 
(Iliadis et al. 2010; § 2.2) and for the temper- 
atures encountered in convective core H-burning, 
all of these are based on experimental data. The 
uncertainties quoted are at the 1-cr level for lognor- 
mal probability density functions. Judging from 
the entries in Tab. [6l none of these reactions 
will impact the abundance of ^^Al^ if their rates 
are varied within a factor of 2 from their rec- 
ommended values. Given the quoted uncertain- 
ties, a factor of 2 corresponds to confidence in- 
tervals of 97.9% for 25Mg(p,7)26A19, 98.4% for 
25Mg(p,7)26Al", 99.3% for ^^AW{p,-f)^'^Si and « 
100% for i'^0(p,7)i'^F. In other words it is unlikely 
that the rates for any of these reactions will be 
a factor of 2 away from the recommended val- 
ues. This points to the utility of uncertainties 
with statistical significance. Therefore, we con- 
clude that the rates for the reactions that deter- 
mine the abundance of ^^Al^ during convective 



core H-burning are known with sufficient preci- 
sion. 

4. Summary 

We presented a comprehensive investigation of 
the impact of nuclear reaction rate uncertainties 
on 2^A1 production in massive stars. In such stars, 
2^A1 is likely produced in three distinct sites: (i) 
during core collapse via explosive Ne/C burning; 
(ii) during pre-supernova stages in the C/Ne con- 
vective shell, where a fraction of the ^^Al sur- 
vives the subsequent explosion and is ejected into 
the interstellar medium; and (iii) in Wolf-Rayet 
stars, that experience such a strong mass loss that 
even layers located within the H convective core, 
hence significantly enriched in ^^Al, are ejected 
into the interstellar medium. These ^^Al produc- 
tion mechanisms were recently analyzed in detail 
by Limongi & Chieffi (2006). From their stel- 
lar evolution models, we extracted representative 
temperature-density-time profiles and executed a 
large number of post-processing reaction network 
sensitivity calculations. The general strategy con- 
sisted of varying the rates of many reactions in- 
dividually by different factors (in this work, 10, 
2, 0.5 and 0.1) and to analyze the impact of each 
individual reaction rate change on the final ^^Al 
yields. Our results are important for quantifying 
the infiuence of current reaction rate uncertainties 
on predicted ^^Al yields, and for the motivation of 
future laboratory measurements. 

There are a number of novel aspects about 
the present work. First, we employed a new- 
generation library of nuclear reaction and weak in- 
teraction rates, called STARLIB. This library con- 
tains a recent evaluation of experimental Monte 
Carlo reaction rates (Iliadis et al. 2010). Besides 
recommended reaction rates for a grid of temper- 
ature values between 1 MK and 10 GK, the li- 
brary includes in addition for many reactions the 
rate uncertainty factor at each temperature. This 
work represents the first application of STAR- 
LIB. Second, we carefully investigate the equili- 
bration effects of 2^A1. At least two species of Al 
take part in the nucleosynthesis, the ground state 
and the isomeric state. In all previous massive 
star investigations, either a single species or two 
species of ^^Al were taken into account, depend- 
ing on whether thermal equilibrium is achieved or 
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Table 6 

Factor changes of final ^^Al^ abundance resulting from reaction rate variations for 

CONVECTIVE CORE H BURNING"" , ASSUMING TWO SPECIES OF ^^Al 



Reaction'' 




Rate multiplied by 




Source*^ 


Uncertainty*^ 


100 


10 2 


0.5 0.1 


0.01 


26Alff(p,7)27Si 


0.0035 


0.55 ••• 






1110 


31% 


25Mg(p,7)26A19 


1.20 


1.20 


• • • 0.33 


0.039 


1110 


35% 


25Mg(p,7)26Al'" 


0.049 


0.37 ••• 


••• 1.20 


1.20 


1110 


35% 


i60(p,7)"F 








1.70 


1110 


7% 



^The temperature-density-time profile is extracted from a stellar evolution calculation 
of a 8OM0 star with initial solar metallicity, see Limongi & ChiefR (2006). 

''In total, the rates of 26 different reactions were varied. Listed are only those reactions 
whose rate changes have the strongest effect on the ^^AF yield. All other rate changes, as 
well as those labeled by produced abundance changes of less than 20%. The reactions 
are listed in approximate order of importance. 

■^Reaction rate reference: (illO) Iliadis et al. (2010). 

■^Reaction rate uncertainty near a temperature of 0.09 GK, at the end of the calculation. 
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Fig. 4. — Abundance evolution (by mass) of ^^Al* 
and ^^Al during explosive Ne/C burning. The T-p 
profile for this post-processing network simulation 
is shown in Fig. The calculation assumes a 
single species of (thermalized) ^^Al. 



not. These are two extreme assumptions and in a 
hot stellar plasma the ground and isomeric state 
may "communicate" via 7-ray transitions involv- 
ing higher-lying ^^Al levels. 

Some of our results are summarized in Tab. 
[71 listing those nuclear reactions that signifi- 
cantly impact ^^Al synthesis in massive stars. 
The reactions are listed in approximate or- 
der of importance. The reader should consult 
Tabs. [21 m and [H for detailed results. Par- 
ticularly the first five reactions, ^^Al(n,p)^^Mg, 
25Mg(a,n)28Si, 24Mg(n,7)25Mg, 23Na(a,p)26Mg 
and 2^Al(n,a)2^Na, should be prime targets for 
future measurements. The approximate temper- 
ature range near which the rate needs to be im- 
proved (w 2.3 GK for explosive Ne/C burning, 

1.4 GK for convective shell C/Ne burning), as 
well as the current literature source of a particular 
rate, is also given in the table. For those five reac- 
tions we argued in ii l3.2.3| 13. 3. 31 and Appendix [Q 
that the current rate uncertainties at astrophysi- 
cally important temperatures amount to about a 
factor of 2. The sensitivity of ^^Al production to 
rate variations of these reactions can be estimated 
from Tabs. [2 and [H a factor of 2 variation in 
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Fig. 5. — Abundance evolution (by mass) of ^^Al 
during explosive Ne/C burning. The T-p pro- 
file for this post-processing network simulation is 
shown in Fig. [2] The calculation assumes five 
species of ^^Al: ground state {g), isomeric state 
(m), and three excited levels (a, 6, c); see Fig. 
[TJ The communication of the different ^^Al lev- 
els via 7-ray transitions is explicitly taken into 
account. (Top) Abundance evolution of differ- 
ent ^^Al species; (Bottom) Numerically simulated 
fraction of total ^^Al abundance that resides in the 
isomeric state. The curve is indistinguishable from 
the one calculated assuming a Boltzmann distribu- 
tion (i.e., thermal equilibrium). 
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Fig. 6. — Temperature-density evolution for con- 
vective shell C/Ne burning. The results were ob- 
tained from a model of a 6QMq star (Limongi & 
ChiefR 2006), but the time scale is shortened in the 
present work (see text). The profile approximates 
the evolution of the hottest and deepest zone of the 
C/Ne convective shell. For comparison, the circles 
indicate the temperature and density values that 
are directly obtained from the stellar evolution cal- 
culations. Time increases from left to right. 
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their rates changes the final ^^Al mass fraction by 
factors of 1.7, 1.9, 1.6, 1.3 and 1.3, respectively. 
Thus we conclude that the uncertainty of the ^^Al 
yield predicted by the massive star models explored 
here amounts to about a factor of 3. This result is 
obtained on the basis of nuclear physics uncertain- 
ties alone and should be considered together with 
other uncertainties inherent in the stellar models, 
such as mixing, mass loss and rotation. 

We do not list any reactions for core H burn- 
ing in Tab. |6l mainly because in this case reac- 
tion rate variations have only a small effect on the 
^^Al yield. Here, the most important reaction is 
^^Al(p,7)^^Si, but even a factor of 10 change in 
this rate near « 90 MK has only a modest impact 
on the ^^Al yield (see Tab. [6]). Of course, new ex- 
perimental results for ^^Al(p,7)^^Si are useful in 
any case. 

We carefully examined the issue of ^^Al equili- 
bration for each of the three nucleosynthesis sites 
mentioned above. Two series of post-processing 
calculations were performed and the resulting ^^Al 
yields were compared: one assuming either a sin- 
gle or two separate ^^Al species, depending on the 
temperature regime, and one where the commu- 
nication between ground and isomeric states was 
explicitly taken into account. For the latter case, 
no artificial assumptions about the equilibration 
of ^^Al are made, but additional ^^Al species (i.e., 
levels at 417, 1058 and 2070 keV; see Fig. ^ 
were taken into account in the reaction network. 
We found that the equilibration of ^^Al levels in 
any of the massive star sites investigated here has 
only minor effects on the ■^^Al yields. The rea- 
son is that in explosive Ne/C burning and con- 
vective shell C/Ne burning the temperatures are 
sufficiently high to ensure thermal equilibration of 
^^Al, while in core H burning the temperatures are 
never high enough to facilitate communication of 
the ground and isomeric state via thermal excita- 
tions. We also verified that current uncertainties 
in some unmeasured ^^Al 7-ray transition rates do 
not significantly impact the predicted nucleosyn- 
thesis yields. 

For the interested reader we provide detailed 
comments on the status of certain reactions, in- 
cluding i2c(i2c,n)23Mg, 23Na(a,p)26Mg, 25Mg(a,n)28Si, 
26Al"'(p,7)27Si, 26Al(n,p)26Mg and 26Al(n,a)23Na. 
For the latter two, particularly important, reac- 
tions we provide new rate estimates, which will be 



presented in more detail in a forthcoming publi- 
cation (Oginni et al., in print). 

This work was supported in part by the U.S. 
Department of Energy under Contract No. DE- 
FG02-97ER41041. 
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Table 7 

Summary of nuclear reactions that impact ^^Al production in massive stars'' , assuming 

THERMAL EQUILIBRIUM FOR ^^Al 



Reaction 


Sitei^ 


Temperature'^ 


Source'* 


26Al*(n,p)26Mg 


xNe/C; C/Ne 


« 2.3; « 1.4 


present 


25Mg(a,n)28Si 


xNe/C 


w 2.3 


nacr 


24Mg(n,7)25Mg 


xNe/C; C/Ne 


w 2.3; « 1.4 


ka02 


23Na(a,p)26Mg 


C/Ne 


« 1.4 


rath 


26Al*(n,a)23Na 


xNe/C; C/Ne 


« 2.3; w 1.4 


present 


27Al(a,p)30Si 


xNe/C 


« 2.3 


rath 


29Si(a,n)32s 


xNe/C 


« 2.3 


rath 


26Mg(a,n)29Si 


C/Ne 


« 1.4 


nacr 



^In approximate order of importance; for fuU resuhs, see 
Tabs. H HI and m 

^Site of 2^A1 synthesis in massive star; the labels "xNe/C" 
and "C/Ne" refer to explosive Ne/C burning and convective 
shell C/Ne burning, respectively. 

'^Temperature (in units of GK) near which most of ^^Al 
production occurs in given site. 

'^For reference labels, see Tabs. [2]and|4l 
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Fig. 8. — Abundance evolution (by mass) of 
and ^^Al during convective shell C/Ne burni: 
The T-p profile for this post-processing netwi 
simulation is shown in Fig. [5] The calculation as- 
sumes a single species of (thermalized) ^^Al. Time 
increases from left to right. 



: Convective shell 




C/Ne burning 






a 
y 




a 


r ^ IP* 


m 


, f 






b 




c i 


0.15 0.14 0.13 0.12 


0.11 



Mass fraction ^^C 



10' 




< 



-\Q-3 I , , , , , 1 

0.15 0.14 0.13 0.12 0.11 0.1 

Mass fraction ^^C 

Fig. 9. — Abundance evolution (by mass) of ^^Al 
during convective shell C/Ne burning. The T-p 
profile for this post-processing network simulation 
is shown in Fig. [S] The calculation assumes five 
species of ^^Al: ground state (g), isomeric state 
(m), and three excited levels (a, 6, c); see Fig. 
[TJ The communication of the different ^^Al lev- 
els via 7-ray transitions is explicitly taken into 
account. (Top) Abundance evolution of differ- 
ent ^^Al species; (Bottom) Numerically simulated 
fraction of total ^^Al abundance that resides in the 
isomeric state. The curve is indistinguishable from 
the one calculated assuming a Boltzmann distribu- 
tion (i.e., thermal equilibrium). 
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Fig. 10. — Temperature-density evolution for con- 
vective core H burning. The results "were obtained 
from a model of a 80Mq star (Limongi & Chi- 
effi 2006), but the time scale is shortened in the 
present work (see text). The circles indicate the 
temperature and density values that are directly 
obtained from the above stellar evolution calcula- 
tions. Time increases from left to right. 
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Fig. 12. — Abundance evolution (by mass) of 
^^Al^ and ^''Al during convective core H burning. 
The T-p profile for this post-processing simulation 
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A. /3- AND 7-DECAY RATES OF ^^Al, LEVELS 



The decay constants for j3- and 7-decay of ^^Al levels, in units of s are listed in Tab. [8] for the 
temperature range of relevance in the present work. The labels ^^AF, ^^Al™, ^^Al°, ^^Al^ and ^^AF refer 
to the levels aX ^ Q keV (5+; ground state), 228 keV (0+; isomeric state), 417 keV (3+), 1058 keV (1+) 
and 2070 keV (2+), respectively (see Fig. [1]). 

The entries in the second column refer to the /3-decay of ^^Al* (ground and isomeric state in thermal 
equilibrium) to the daughter ^^Mg. The decay constant is only listed for temperatures above T — 0.4 
GK since for lower temperatures thermal equilibrium is not achieved. The values are calculated from 
\{^^k\*^'^^Mg) = 9.93 X 10-3 e-2-646/T9 g-i^ ^heve Tg is the temperature in GK (see Ward & Fowler 
1980, Iliadis 2007). This expression, which takes only the ground and isomeric state into account, is valid 
for temperatures and densities below 5 GK and 10^ g/cm'^, respectively. In the temperature and density 
regimes considered here, the results are in good agreement with the more extensive calculations of Oda et 
al. (1994). Note that the original REACLIB fit of this particular decay rate is off by «20-80% at T = 1 - 3 
GK. 

The following ^-decay constants are not listed in the table since they are constant for the temperature 
grid shown here: 

A(26A19 ^26]y[g)^3 069x10-14 g^l 

A(26A1™ ^26^g)=l. 092x10-1 
A(26AP ^26i^g)^j^ 283x10-4 s-i 

xl^Gj^^b _^>26^g)^9 630x10-2 g-i 

The first two values are computed from measured laboratory half-lifes (Audi et al. 2003), while the latter 
two values are obtained from shell model calculations (Kajino et al. 1988). 

The following 7-ray decay constants are not listed in the table since they are nearly constant for the 
temperature grid shown here: 
X{^^Af ^■26A1™)=2. 780x1013 g-^ 
A(26A1= ->26A1™)=1. 500x1012 s-i 
A(26A1^ ^.26AP)=7.240xl08 s-i 

A(26AF Ala) = l 070X1013 fi-l 

A(26A1= _^26Aib)=3.770xl0i3 5-1 

All 7-ray decay constants given above and listed in Tab. |8] are calculated from experimental information 
(Endt 1990), except those connecting the levels 26^1™ o^'^Al" and 26^1" .(^.26ai&^ which have been found 
from shell model calculations (Runkle, Champagne & Engel 2001; see also Fig. [1]). 
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B. REACTION RATES INVOLVING ^saL*, ^saL^ and ^saL 



Information on the reaction rates involving the production and destruction of ^^Al species in our network 
is given in Tab. [SI The second row hsts the source of the rates: "illO" (experimental Monte Carlo rates from 
Iliadis et al. 2010); "rath" (theoretical Hauser-Feshbach rates from Rauscher & Thielemann 2000); "present" 
(hybrid rate, see below) . We argued in § 12.31 that it is important to ensure internal consistency of the rates 
used. For example, for the first three listed reactions (^^Mg+p) the rates are based on the same nuclear 
physics input and are thus consistent. Similar arguments apply to the following two reactions (^^AP+p). 

However, the situation for the ^^Al™(p,7)^^Si reaction is a diflferent matter. Rates have been estimated 
in Caughlan & Fowler (1988) and in Angulo et al. (1999), while initial experimental studies are reported in 
Deibel et al. (2009) and Lotay et al. (2009). There can be no doubt that this rate is highly uncertain at 
present (see Appendix IC. 41) . In the absence of a better procedure, we approximated the ^^Al™(p,7)^^Si rate 
by the (experimental) ground state rate (column 6 in Tab. [3]). Note that these two rates were predicted 
by Caughlan & Fowler (1988) to be similar within a factor of «5. Our assumption should be regarded as a 
starting point for exploring the effects of ^^Al™(p,7)^^Si reaction rate variations. 

For the ^^Al*(n,a)^^Na reaction we use a hybrid rate, which is based on experimental information from 
De Smet et al. (2007) at T < 0.1 GK, and on Hauser-Feshbach results from Rauscher & Thielemann (2000) 
at higher temperatures (see Appendix I C. 5 ( the rate is listed in the last column of Tab. |9]). The predicted 
stellar enhancement factors are relatively small (43% at 2.5 GK) and, therefore, we adopt these rates also 
for the ^^AP(n,a)^'^Na reaction. Furthermore, in order to ensure internal consistency, we approximated 
the ^^Al™(n,a)^^Na rate by the thermalized rate. Note that the thermalized and isomeric state rates were 
predicted by Caughlan & Fowler (1988) to be similar within a factor of !=a7. Again, our assumptions serve 
as starting points to explore the effects of ^^AP(n,a)^'^Na reaction rate variations. A similar procedure has 
been followed for the ^®AP(n,p)^^Mg reaction rates (see column 7 of Tab. [9] and Appendix IC.6|) . 

The rates of the following reactions are not explicitly listed in Tab. [SI ^^Na(a,7)^^Al*, ^^Al(n,7)^^Al*, 
23Mg(a,p)26Al*, 26Si(n,p)26Al*, 29p(n,a)26Al*, 26Al*(n,7)27Al, A\\a,-f)^°P and 26Al*(a,p)29Si. For these 
we used the Hauser-Feshbach rates of Rauscher & Thielemann (2000), except the rates of ^^Al*(n,7)^''Al, 
which were adopted from the KADoNiS vO.2 evaluation (Dillmann et al. 2006). Since the stellar enhancement 
factors are predicted to be small, we also adopted these results for the respective rates involving ^^Al^. Note 
that the corresponding reactions involving 26 Al'" are absent in the original REACLIB. We disregarded these 
as well, on the grounds that their net abundance flows (for ^^Al*) in our network calculations are at least 3 
orders of magnitude smaller than the maximum flow. 

For all forward reactions discussed above, the corresponding reverse reaction rates are also implemented 
in our library. Note that in these cases we estimated the reverse rates using the principle of detailed balance 
assuming a single (ground or isomeric) level in ^^Al only. The proper procedure would have been to apply 
detailed balance to the forward rate involving the quasi-equilibrium cluster of levels in thermal equilibrium 
with the ^^Al level in question. We are not aware that this information has been given anywhere before and 
believe that the effects of our approximation are very small. 
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Table 8 

Decay constants of ^^Al levels (in s^^) 



T (GK) 






26^jo _j.26^[g 


26A1™ ^26 


26 Al" _i.26Al'" 


26 Al™ ^2'=A1'' 


26A1™ ->26ai'^ 


26Ar -^28^j6 


26 Ar ^2«A1<^ 


26^1^ ^26ai<: 


0.04 




1.12E-44 


5.56E+08 


7.11E-25 


6.24E-02 












0.05 




3.43E-34 


5.56E+08 


4.07E-20 


6.24E-02 












0.06 




3.45E-27 


5.56E+08 


6.05E-17 


6.24E-02 












0.07 




3.46E-22 


5.56E+08 


l.llE-14 


6.24E-02 






2.21E-38 






0.08 




1.95E-18 


5.56E+08 


5.58E-13 


6.24E-02 


4.45E-39 




1.30E-32 






0.09 




1.62E-15 


5.56E+08 


1.17E-11 


6.24E-02 


2.86E-33 




3.97E-28 




1.37E-43 


0.10 




3.49E-13 


5.56E+08 


1.33E-10 


6.24E-02 


1.27E-28 




1.54E-24 




6.36E-38 


0.11 




2.83E-11 


5.56E+08 


9.77E-10 


6.24E-02 


8.03E-25 




1.33E-21 




2.75E-33 


0.12 




l.lOE-09 


5.56E+08 


5.14E-09 


6.24E-02 


1.18E-21 




3.73E-19 




2.01E-29 


0.13 




2.45E-08 


5.56E+08 


2.09E-08 


6.24E-02 


5.67E-19 




4.40E-17 




3.72E-26 


0.14 




3.50E-07 


5.56E+08 


6.98E-08 


6.24E-02 


1.13E-16 




2.62E-15 




2.36E-23 


0.15 




3.50E-06 


5.56E+08 


1.98E-07 


6.24E-02 


l.lOE-14 




9.04E-14 


2.24E-43 


6.33E-21 


0.16 




2.63E-05 


5.56E+08 


4.94E-07 


6.24E-02 


6.09E-13 


1.40E-45 


2.00E-12 


6.64E-40 


8.44E-19 


0.18 




7.56E-04 


5.56E+08 


2.26E-06 


6.24E-02 


4.89E-10 


2.05E-39 


3.51E-10 


4.04E-34 


2.93E-15 


0.20 




l.lOE-02 


5.56E+08 


7.63E-06 


6.24E-02 


1.03E-07 


2.94E-34 


2.19E-08 


1.71E-29 


2.00E-12 


0.25 




1.40E+00 


5.56E+08 


6.83E-05 


6.24E-02 


1.56E-03 


5.62E-25 


3.72E-05 


3.66E-21 


2.51E-07 


0.30 




3.52E+01 


5.56E+08 


2.94E-04 


6.24E-02 


9.59E-01 


8.66E-19 


5.29E-03 


1.31E-15 


6.31E-04 


0.35 




3.52E+02 


5.56E+08 


8.37E-04 


6.25E-02 


9.40E+01 


2.28E-14 


1.83E-01 


1.21E-11 


1.69E-01 


0.40 




1.9oE+U3 


5.5DE+U8 


1.8oE-U3 


D.27E-U2 


2.y3E+Uo 


4.7UE-il 


z.dUE+UU 


1.15E-08 


1.12E-(-Ui 


0.45 


2.78E-05 


7.60E+03 


5.56E+08 


3.38E-03 


6.29E-02 


4.25E+04 


1.78E-08 


2.06E+01 


2.36E-06 


2.92E+02 


0.50 


5.00E-05 


2.23E+04 


5.56E+08 


5.53E-03 


6.32E-02 


3.61E+05 


2.05E-06 


1.07E+02 


1.67E-04 


3.97E+03 


0.60 


1.21E-04 


1.12E+05 


5.56E+08 


1.20E-02 


6.41E-02 


8.94E+06 


2.55E-03 


1.28E+03 


l.OOE-01 


1.99E+05 


0.70 


2.27E-04 


3.54E+05 


5.57E+08 


2.00E-02 


6.53E-02 


8.85E+07 


4.13E-01 


7.53E+03 


9.63E+00 


3.26E+06 


0.80 


3.64E-04 


8.40E+05 


5.57E+08 


3.00E-02 


6.67E-02 


4.94E+08 


1.88E+01 


2.84E+04 


2.96E+02 


2.65E+07 


0.90 


5.25E-04 


1.65E+06 


5.59E+08 


4.20E-02 


6.84E-02 


1.88E+09 


3.65E+02 


7.98E+04 


4.24E+03 


1.35E+08 


1.00 


7.04E-04 


2.83E+06 


5.60E+08 


5.50E-02 


7.02E-02 


5.49E+09 


3.92E+03 


1.82E+05 


3.58E+04 


4.99E+08 


1.25 


1.20E-03 


7.54E+06 


5.68E+08 


9.20E-02 


7.55E-02 


3.77E+10 


2.82E+05 


8.10E+05 


1.66E+06 


5.23E+09 


1.50 


1.70E-03 


1.47E+07 


5.79E+08 


1.32E-01 


8.13E-02 


1.36E+11 


4.87E+06 


2.19E+06 


2.14E+07 


2.50E+10 


1.75 


2.19E-03 


2.38E+07 


5.93E+08 


1.75E-01 


8.74E-02 


3.42E+11 


3.73E+07 


4.49E+06 


1.33E+08 


7.66E+10 


2.00 


2.64E-03 


3.46E+07 


6.10E+08 


2.20E-01 


9.38E-02 


6.82E+11 


1.72E+08 


7.71E+06 


5.23E+08 


1.78E+11 


2.50 


3.45E-03 


5.98E+07 


6.50E+08 


3.12E-01 


1.07E-01 


1.81E+12 


1.45E+09 


1.67E+07 


3.56E+09 


5.78E+11 



Note. — Tabic \s\ is published in the electronic edition of the Astrophysical Journal. 



Table 9 

Stellar reaction rates''' Na{ijv) involving ^^Al (in cm^ mol~^ s"-^) 



to 



T (GK) 


=^Mg(p,7)26Al* 


2^Mg(p,7)^=Alf 


==Mg(p,7)^=Al"' 


=«'Al*(p,7)=^Si 


=«'Alf(p,7)=''Si 


=«Al*(n,p)=6Mg 


=«Al*(n,a)=3Na 




illO'' 


illO'' 


illO'' 


illO'' 


illO''''^ 


present'^'''' 


present^'-^ 


0.015 


1.07O-24 


8.71e-25 


2.04O-25 


4.07C-31 


4.07e-31 


2.98O+06 


4.39O+06 


0.016 


1.560-23 


1.270-23 


2.960-24 


4.57O-30 


4.57e-30 


3.4O0+O6 


5.2O0+O6 


0.018 


1.350-21 


1.09O-21 


2.550-22 


6.05O-28 


6.05e-28 


4.19O+06 


6.77O+06 


0.020 


4.69e-20 


3.81e-20 


8.89e-21 


4.78e-26 


4.78e-26 


4.9O0+O6 


8.17O-I-06 


0.025 


2.67e-17 


2.16e-17 


5.05e-18 


2.73e-22 


2.73e-22 


6.37e4-06 


1.09e-|-07 


0.030 


1.74e-15 


1.41e-15 


3.30e-16 


1.49e-19 


1.49e-19 


7.53e-|-06 


1.25e-|-07 


0.040 


3. 000-13 


2.440-13 


5.68e-14 


7.35e-16 


7.35e-16 


9.41e+06 


1.39O-F07 


0.050 


6.820-12 


5.570-12 


1.260-12 


1.25C-13 


1.250-13 


I.II0+07 


I.4I0+O7 


0.060 


6.560-11 


5.430-11 


1.150-11 


3.91e-12 


3.910-12 


1.27O+07 


1.39O+07 


0.070 


3.99O-10 


3.33O-10 


6.670-11 


5.71e-ll 


5.710-11 


1.42O+07 


1.37O+07 


0.080 


1.72e-09 


1.45e-09 


2.81e-10 


7.80e-10 


7.80e-10 


1.58e+07 


1.34O-I-07 


0.090 


5.73e-09 


4.85e-09 


9.35e-10 


9.54e-09 


9.54e-09 


1.74e-|-07 


1.31e-|-07 


0.100 


1.60e-08 


1.34e-08 


2.72e-09 


8.17e-08 


8.17e-08 


1.91e-|-07 


1.28e-|-07 


0.110 


4.14e-08 


3.40e-08 


7.78e-09 


4.94e-07 


4.94e-07 


2.07e-F07 


1.26e-|-07 


0.120 


1.17e-07 


9.36e-08 


2.38e-08 


2.24e-06 


2.24e-06 


2.24e-F07 


1.24e-|-07 


0.130 


4.05e-07 


3.25e-07 


8.12e-08 


8.10e-06 


8.IO0-O6 


2.41e+07 


1.23O+07 


0.140 


1.62e-06 


1.32e-06 


2.93e-07 


2.45e-05 


2.45O-05 


2.58O+07 


I.2I0+O7 


0.150 


6.40O-06 


5.34e-06 


1.05O-06 


6.41e-05 


6.4I0-O5 


2.74e+07 


I.2O0+O7 


0.160 


2.31O-05 


1.95e-05 


3.53O-06 


1.50e-04 


1. 500-04 


2.91e+07 


1.19O+07 


0.180 


2.10e-04 


1.79e-04 


3.02e-05 


6.38e-04 


6.38e-04 


3.24e-t-07 


1.18e-|-07 


0.200 


1.26e-03 


1.08e-03 


1.79e-04 


2.13e-03 


2.13e-03 


3.57e4-07 


1.17e-|-07 


0.250 


3.20e-02 


2.73e-02 


4.68e-03 


2.30e-02 


2.30e-02 


4.36e-F07 


1.17e-|-07 


0.300 


2.76e-01 


2.33e-01 


4.25e-02 


1.40e-01 


1.40e-01 


5.10e-F07 


1.17e-|-07 


0.350 


1.29e-l-00 


1.08e-|-00 


2.10e-01 


5.74e-01 


5.74e-01 


5.77e-|-07 


1.19e-|-07 


0.400 


4. iDC-|-UO 


o.45e+00 


7.06C-01 


i.74e-t-U0 


1 ^ A ^ 1 nf\ 


6.38e-|-07 


1.21o-t-07 


0.450 


1.04C+01 


8.55O+00 


1.85c-t-00 


4.19e-|-00 


4.19o-f-00 


6.93O+07 


1.23O-I-07 


0.500 


2.19C-I-01 


1.78O+01 


4.06O-I-00 


8.50C-I-00 


8.50e-|-00 


7.43O+07 


1.26O-I-07 


0.600 


6.84O-I-01 


5.46O+01 


1.37e-|-01 


2.45C-I-01 


2.45e-|-01 


8.27O+07 


1.32O-I-07 


0.700 


1.58e-l-02 


1.24e-|-02 


3.38e-|-01 


5.14e-l-01 


5.14e-|-01 


8.95e+07 


1.39e-|-07 


0.800 


2.99e-|-02 


2.31e4-02 


6.80e-|-01 


8.96e-|-01 


8.87e-|-01 


9.51e4-07 


1.46e-|-07 


0.900 


4.98e-|-02 


3.78e-F02 


1.19e-|-02 


1.36e-|-02 


1.35e-|-02 


9.98e-F07 


1.53e-|-07 


1.000 


7.53e-|-02 


5.65e-F02 


1.88e-|-02 


1.92e-|-02 


1.88e-|-02 


1.04e-F08 


1.61e-|-07 


1.250 


1.62e-|-03 


1.18e-F03 


4.38e-|-02 


3.59e-|-02 


3.45e-l-02 


1.12e-F08 


1.82e-|-07 


1.500 


2.75e-|-03 


1.98e+03 


7.88e-|-02 


5.57e-|-02 


5.25c-(-02 


1.19e+08 


2.04e-|-07 


1.750 


4.09e-|-03 


2.89O+03 


1.22o-(-03 


7.84o-f02 


7.19O+02 


1.25e+08 


2.29O-I-07 


2.000 


5.56e-|-03 


3.90e-|-03 


1.71e-|-03 


1.14e-|-03 


l.Ole-l-03 


1.31e-|-08 


2.55e-|-07 


2.500 


8.66e-|-03 


6.08e-F03 


2.79e-|-03 


1.99e-|-03 


1.67e-|-03 


1.46e-F08 


3.14e-|-07 



^The rates for only some reactions are listed; other reactions are discussed in the text. All rates given here (except column 6) account for thermal 
target excitations. 

'^Experimental Monte Carlo rates of Iliadis et al. (2010). 

•^Same rate is used for ^^Al^ (p,7)2^Si. 

"^Hybrid rate: at T < 0.2 GK from experiment of Koehler et al. (1997); at T > 0.2 GK from Hauser-Feshbach model of Rauscher & Thielemann 
(2000). 



'^Hybrid rate: at T < 0.1 GK from experiment of De Smet ct al. (2007); at T > 0.1 GK from Hauser-Feshbach model of Rauschcr & Thielemann 
(2000). 

*Samc rate is used for '^^Al^ and ^"^Al"^. 

Note. — Table[9]is published in the electronic edition of the Astrophysical Journal. 
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C. DISCUSSION OF SPECIFIC REACTION RATES 



The references of the reaction rates used in the present work are provided in the tables above and the 
reader is referred to these sources for details. Decays and reactions involving levels of ^^Al are discussed in 
Appendices |^ and [B] respectively. For some specific reaction rates that are neither evaluated in Iliadis et 
al. (2010) nor in Angulo et al. (1999), we provide more information below. 

C.l. i2c(i2c,n)23Mg (Q = -2.598 MeV) 

We calculated the reaction rate from the total ^^C+^^C rate (i.e., summed over all exit channels) and the 
neutron branching ratio. For T = 1.25 GK, representing the average temperature of convective shell C/Ne 
burning of the 60 Mq model explored in the present work, the Gamow peak extends over a center-of-mass 
energy range from 2.6 MeV (i.e., the threshold energy for this endothermic reaction) to about 3.4 MeV. In 
this energy region the total ^^C+^^C cross section has been measured (Costantini et al. 2009, and references 
therein). For the neutron branching ratio, as a function of temperature in the range of T = 0.5 — 5.0 GK, we 
adopted the values from Dayras, Switkowski and Woosley (1977). In the latter work, the ^^C(^^C,n)^^Mg 
reaction has been measured down to an energy of 3.54 MeV and, with the aid of statistical model calculations, 
the results were extrapolated down to the threshold energy in order to extract the neutron branching ratio. 
The overall reaction rate uncertainty is difficult to quantify at present. Clearly, a new measurement of the 
-'^^C(^^C,n)^^Mg reaction at lower energies is desirable. 

The ^^C('^^C,n)^^Mg reaction rates employed in various rate libraries are inconsistent with each other. 
An impression can be obtained from Fig. 1131 The black line represents the ratio of the rate used in the 
Basel version of the REACLIB (nucastro.org/reaclib. html^reaclib) to the present rate. The hatched area 
marks the temperature range of convective shell C/Ne burning explored in the present work. Surprisingly, 
the ratio amounts to more than an order of magnitude near T = 1.25 GK. The reason may be that the 
^^C(^^C,n)^^Mg rate in the original REACLIB library was erroneously obtained using a neutron branching 
ratio of zero for T < 1.75 GK that is listed in Caughlan & Fowler (1988). The blue line in Fig. [T51 represents 
the ratio of the rate used in the MSU version of the REACLIB (groups.nscl.msu.edu/jina/reaclib/db/) to 
the present rate. The ratio amounts to a factor of 2 near T = 1.25 GK. The reason may perhaps be that 
their fitted rate deviates from the actual rate. For the interested reader, we provide below a rate fit (in 
REACLIB format) that reproduces the present and best currently available rate within 10%: 

-0.527166E+02 
-0.344948E+02 

0.140849E+03 
-0.878184E+02 

0.371377E+01 
-0.338673E+00 

0.736030E+02 

C.2. 23Na(a,p)26Mg (Q = 1.821 MeV) 

A direct measurement of this reaction has been reported in Whitmire & Davids (1974). They bombarded a 
target, fabricated by evaporating NaCl onto a thick Cu backing, with a-particles (Ecm = 2.0 — 3.1 MeV) and 
measured the emitted protons populating the ground and first excited states in ^^Mg. In total, they report 
the strengths of 39 resonances. There are a number of reasons that warrant a re-measurement of this reaction. 
Most importantly, the strengths have been determined relative to an absolute strength measurement for the 
Ejaf) = 3051 keV resonance, assuming that the stoichiometry of their NaCl target amounts to 1:1. This issue 
has been discussed in detail by Rowland et al. (2002) in connection with the ^^Na(p,a)^*'Ne reaction, where it 
was shown that during proton bombardment a NaCl target quickly changes its stoichiometry to 5:3, resulting 
in a significant change in thick-target yield (and in the derived resonance strength). This problem is certainly 
aggravated when using an a-particle beam incident on a NaCl target. There were other problems in the 
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analysis of Whitmire & Davids (1974): (i) protons populating higher-lying final states could only be resolved 
at the highest measured energies; (ii) the uncertainty in the resonance energies is relatively large, amounting 
to «10 keV; (iii) the assumed uncertainties of 1% for the stopping cross sections of a-particles in Na and CI 
seem unreasonably small. Furthermore, for temperatures typical of convective shell C/Ne burning {T « 1.25 
GK) the Gamow peak is covering a center-of-mass energy range of Ecm = 1-2 — 2.2 MeV, i.e., significantly 
below the energy range covered by experiment. Considering the substantial uncertainties involved, we prefer 
to use for the reaction rate the estimate provided by the Hauser-Feshbach model (Rauscher & Thielemann 
2000) . Clearly, an improved measurement of this reaction is called for. 

C.3. 25]v[g(Q,^n)28si (Q ^ 2.654 MeV) 

Direct measurements of the ^^Mg(a,n)^^Si and ^^Mg(a,n7)^^Si reactions have been reported in Van der 
Zwan and Geiger (1981), Anderson et al. (1983) and Wieland (1995). The reaction rate recommended 
by the NACRE collaboration (Angulo et al. 1999) is exclusively based on these data sets. The NACRE 
collaboration also reports the experimental S-factors for the three references. However, below an energy of 3 
MeV the data from Van der Zwan and Geiger (1981) and from Anderson et al. (1983) have been disregarded 
by NACRE. As a result, in this energy range, their rate is based exclusively on the unpublished work of 
Wieland (1995), since it is argued in Angulo et al. (1999) that at lower energy background contributions in 
the earlier works dominate the neutron yield. This conclusion is only partially correct since, for example, 
Anderson et al. (1983) have also measured the ^^Mg(a,n7)^^Si reaction, which shows much less background 
compared to the {a,n) reaction. 

The current situation is shown in Fig. [141 The data from Van der Zwan and Geiger (1981) and from 
Anderson et al. (1983) have been extracted as cross sections from the original figures and converted to 
astrophysical S-factors. The data from Wieland (1995) are adopted from Angulo et al. (1999). It is 
interesting to note that the data of Anderson et al. (1983) as reported by the NACRE collaboration disagree 
with the corresponding curve shown in Fig. [14] by a factor of sa2. The reason is presumably that NACRE 
extracted the (a,n) data (which represent only upper limits, as explicitly stated in Anderson et al. 1983) 
instead of the (a,n7) data (which are much less susceptible to background). We conclude from Fig. [T4lthat 
the available data are in reasonable agreement. 

It is also inexplicable why already "...above Tg = 2, H[auser]F[eshbach] rates are used..." in the NACRE 
evaluation. The Gamow peak region at T = 2.3 GK is shown in Fig. [T4| as a hatched horizontal bar. 
Clearly, the Gamow peak region at this temperature, in fact, up to T = 4 GK, is entirely covered by 
experimental data. Thus, there is no reason to use Hauser-Feshbach rates and thereby introduce another 
source of uncertainty. We feel that a proper re-analysis of the existing cross section data will not only improve 
the rate of the ^^Mg(Q!,n)^*Si reaction, but will likely provide better estimates of the rate uncertainty, at 
least up to temperatures of 4 GK. Such an analysis is left for future work. We would also like to reiterate 
that a new measurement of this reaction would be important. 

C.4. 26Al™(p,-y)27si (Q ^ 7.691 MeV) 

Rates for this reaction are listed in Caughlan & Fowler (1988), but it is not apparent from their work how 
the results have been obtained. Presumably their rates were estimated using statistical model calculations 
(see comments in Ward & Fowler 1980). In Angulo et al. (1999), the rate for this reaction was obtained 
by multiplying the experimental (ground state) rate for ^^Als(p,7)^^Si by the ratio of isomeric and ground 
state rates, Na {av)^/NA {<^v)g. The latter ratio was obtained from the Hauser-Feshbach model. However, 
it is clear from our comments in § 13.2.31 that the latter reaction model may not be applicable to ^^Al-|-p. 
Furthermore, the rate in Angulo et al. (1999) is only listed at temperatures of T = 0.018 — 0.4 GK, but the 
rate needs to be known at higher temperatures as well in order to study the equilibration of ^^Al levels (§ 
[231). 

Recently, some new experimental information has been reported by Deibel et al. (2009) and Lotay et al. 
(2009). In the former work, levels in the ^''Si compound nucleus near the proton threshold were populated 
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Fig. 13. — Reaction rate ratio for ^^C(^^C,n)^'^Mg: (black line) ratio of "Basel" rate to present rate; (blue 
line) ratio of "MSU" rate to present rate. The hatched region marks the temperature range of convective 
shell C/Ne burning explored in the present work. 




E (MeV) 

cm 

Fig. 14. — Experimental astrophysical S-factors for ^^Mg(a,n)^^Si: (red) Van der Zwan and Geiger (1981); 
(blue) Anderson et al. (1983); (black) Wieland (1995). For the first two references, the data were extracted 
from the published figures and converted from cross sections to S-factors. For the latter reference, the S- 
factors are adopted from the NACRE evaluation (Angulo et al. 1999). Uncertainty bars have been omitted 
for reasons of clarity. The hatched horizontal bar marks the region of the Gamow peak near T — 2.3 GK, 
the peak temperature achieved in explosive Ne/C burning. 
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in (^He,t) and (^He,a) reaction studies and the subsequent proton decay to the isomeric state was observed 
in coincidence, providing values for excitation energies and proton branching ratios. In the latter work, the 
12(--;^16q reaction was used to measure 7-ray transitions in ^^Si, allowing for a determination of excitation 
energies, J'^-values and level lifetimes. Nevertheless, too much experimental information is still lacking (i.e., 
missing levels, spectroscopic factors, proton partial widths, and resonance strengths) in order to estimate 
this rate reliably over the temperature range of interest. More measurements are clearly in order. 

In the absence of a more reliable estimate, we approximated in this work the ^^Al™(p,7)^^Si rate by the 
(experimental) ground state rate (see comments in Appendix |BJ . Our assumption is a starting point for 
exploring the effects of ^^Al™(p,7)^^Si reaction rate variations. 

C.5. ^^A\(n,a)^^Na (Q ^ 2.966 MeV) 

A direct measurement of the ^^Al(n,Q!o)^'^Na reaction (i.e., for population of the ^^Na ground sate) has 
been reported by Koehler et al. (1997), while De Smet et al. (2007) have measured the ^^Al(n,ao + Q;i)^^Na 
reaction (i.e., for population of the ground and first excited state in ^^Na). The current situation for the 
reaction rates is displayed in Fig. [15] The experimental rate for ^^Al(n,ao)^^Na from Koehler et al. (1997) 
is extracted from their Fig. 4 and is shown as a dashed line. Note that their rates are claimed to be reliable 
only for T < 0.08 GK. The more recent experimental rate of De Smet et al. (2007) is shown as a black solid 
line and was obtained by converting the Maxwellian-averaged cross sections (MACS), shown in their Fig. 8, 
to reaction rates. This rate represents a lower limit above T = 0.26 GK (indicated by the vertical line). The 
theoretical rate, based on the Hauser-Feshbach model, is adopted from Rauscher & Thielemann 2000 and 
is displayed as a solid blue line. A similar theoretical rate has been reported by Goriely, Hilaire & Koning 
(2008). 

We may draw a number of conclusions from the figure. First, the two experimental rates do not agree, 
even if the large rate uncertainty (26%) in the earlier work is taken into account (see discussion in De Smet 
et al. 2007 for the possible source of the discrepancy). Second, taking both the experimental uncertainty of 
the De Smet rate into account, as well as the fact that their rate represents a lower limit near their high- 
temperature cutoff, the agreement with the Hauser-Feshbach rate near T « 0.3 GK is reasonable. Recall 
that for the purposes of the present work the rate is of interest at temperatures between 1.1 GK (convective 
shell C/Ne burning) and 2.3 GK (explosive Ne/C burning), i.e., in a region that has not been covered by 
experiments. Current reaction rate uncertainties are difficult to quantify and new measurements are called 
for. In the absence of more reliable results, we use a hybrid rate consisting of the results from De Smet et 
al. (2007) and from Rauscher & Thielemann (2000) below and above T w 0.1 GK, respectively. The rate is 
listed in the last column of Tab. [9l 

C.6. 26Al(n,p)26Mg (Q = 4.787 MeV) 

The 2^Al(n,p)2^Mg reaction for the transitions to both the ground and first excited ^^Mg state was directly 
measured by Trautvetter et al. (1986) at a number of neutron energies (corresponding to T « 0.36 — 3.6 GK). 
Koehler et al. (1997) measured the 2^Al(n,pi)2^Mg reaction, i.e., for population of the first excited state in 
2^Mg. The latter rates are claimed to be reliable only for T < 0.3 GK. Near the overlap region, T « 0.3 GK, 
the latter rate exceeds the former rate by a factor of 2 (see discussion in Koehler et al. 1997 for the possible 
source of the discrepancy). The disagreement cannot be explained by unaccounted transitions, because the 
Koehler rate exceeds the Trautvetter rate. In any case, it is shown in both Trautvetter et al. (1986) and 
in Skelton, Kavanagh & Sargood (1987) that the 2^Al(n,po)^^Mg reaction rate (i.e., for population of the 
ground state in ^^Mg) is predicted to be much smaller than the (n,pi) rate. 

The current situation for the reaction rates is displayed in Fig. 1161 The experimental rate for 
2^Al(n,pi)2^Mg from the more recent work of Koehler et al. (1997) is extracted from their Fig. 5 and is 
shown as a dashed line. The theoretical 2^Al(n,p)2^Mg rate, based on the Hauser-Feshbach model (Rauscher 
& Thielemann 2000), is displayed as a solid blue line. An almost identical theoretical rate has been reported 
by Goriely, Hilaire & Koning (2008). Considering the experimental uncertainty of the Koehler rate near 
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their high-temperature cutoff (20%), the agreement between experimental and theoretical rates near T w 0.2 
GK seems reasonable (deviation of 40%). Therefore, we adopt a hybrid rate consisting of the results from 
Koehler et al. (1997) and from Rauscher & Thielemann (2000) below and above T sa 0.2 GK, respectively. 
The rate is listed in column 7 of Tab. [5] Note that we prefer the more recent rates from Koehler et al. 
(1997) and Rauscher & Thielemann (2000) over the earlier experimental result of Trautvetter et al. (1986), 
which is displayed as data points in the figure. Our adopted rate exceeds the prediction of Trautvetter et 
al. (1986) by a factor of «3 near T ~ 2.5 GK. It is currently difficult to estimate rate uncertainties and new 
measurements are urgently needed. 
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Fig. 15.— Reaction rates for ^'^A\{n,a)^^'Na,: (black solid line) De Smet et al. (2007); (dashed line) Koehler 
et al. (1997); (blue solid line) Rauscher & Thielemann (2000). The first two rates are based on experimental 
results, while the latter rate is estimated using the Hauser-Feshbach model. Beyond the vertical line, near 
T X 0.26 GK, the experimental rate of De Smet et al. (2007) represents a lower limit. Note that for 
this comparison only, the rates represent "laboratory rates", i.e., they do not account for thermal target 
excitations. 




Fig. 16. — Reaction rates for ^^Al(n,p)^^Mg: (dashed line) Koehler et al. (1997); (blue sohd line) Rauscher & 
Thielemann (2000); squares (Trautvetter et al. 1986). The first (experimental) rate only takes the transition 
to the first excited state in ^^Mg into account, while the third (experimental) rate represents the combined 
transitions to the ground and first excited states in ^^Mg. The second rate is estimated using the Hauser- 
Feshbach theoretical model and includes transitions to all possible final states. Note that for this comparison 
only, the rates represent "laboratory rates", i.e., they do not account for thermal target excitations. 
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